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Summary 


A  unique  characteristics  of  rotorcraft  is  the  dominating  effect  of  the  rotor  down- 
wash  on  the  aerodynamics  of  the  rotor  and  the  airframe.  This  study  addresses  en¬ 
hancements  to  inflow  modeling  that  significantly  improve  the  prediction  of  helicopter 
behavior  across  the  flight  envelope  when  used  in  conjunction  with  a  comprehensive 
rotorcraft  simulation.  The  simulation  includes  interaction  of  a  broad  range  of  phe¬ 
nomena  and  subsystems  including  rotor  dynamics,  rotor  wake,  unsteady  airloads  with 
dynamic  stall,  the  flight  control  system  and  the  propulsion  (including  engine  dynamics 
and  fuel  control). 

Specifically,  the  model  enhancements  that  have  been  made  under  the  Phase  One 
effort  can  be  grouped  into  three  areas. 

•  The  tip  vortex  distortion  effect  due  to  the  hub  rotation  has  been  included  in 
the  finite  state  dynamic  wake  model  formulation  to  overcome  the  poor  predic¬ 
tion  of  rotorcraft  off-axis  response.  Also  a  significant  breakthrough  has  been 
accomplished  in  the  development  of  a  parametric  formulation  of  the  finite  state 
dynamic  wake  model  with  wake  distortion  influence  coefficients  derived  from  a 
distorted  vortex  wake  model  in  both  hover  and  forward  flight. 

•  The  enhanced  formulation  allows  for  better  prediction  of  aerodynamic  interfer¬ 
ence.  This  includes  the  rotor  aerodynamic  interference  on  other  lifting  surfaces 
such  as  the  horizontal  tail,  the  vertical  fin,  and  the  tail  rotor. 

•  A  modification  of  the  finite  state  dynamic  wake  model  has  been  made  to  include 
the  vortex  ring  state.  The  resulting  simulation  is  applicable  across  the  flight 
envelope. 

Investigations  have  been  performed  to  prove  the  feasibility  of  the  approach.  The 
enhanced  model  has  been  validated  against  measured  data  in  hover,  forward  flight 
and  in  the  vortex  ring  state.  The  results  have  shown  a  remarkable  improvement  in 
the  prediction  of  helicopter  response  in  both  hover  and  forward  flight. 
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Chapter  1 
Introduction 


Modern  computer  simulation  and  computer  aided  analysis  technology  have  found 
widespread  applications  in  various  engineering  fields.  Its  potential  for  shortening 
the  rotorcraft  design  cycle,  reducing  cost,  enhancing  flight  test  productivity,  and 
increasing  the  safety  of  flight  operations  is  yet  to  be  explored.  High  fidelity  piloted 
flight  simulation  models  are  required  to  support  land  and  sea-based  flight  testing 
and  operational  analysis.  Such  models  need  to  represent  rotorcraft  accurately  before 
they  can  be  utilized  in  both  engineering  analysis  for  flight  test  support  and  piloted 
flight  simulation  for  training  under  routine  or  emergency  conditions.  Among  other 
things,  a  reliable  prediction  of  rotorcraft  dynamic  response  to  flight  control  inputs 
is  a  fundamental  requirement  in  piloted  flight  simulation.  Although  many  efforts 
have  been  undertaken  to  improve  rotorcraft  mathematical  modeling,  the  capability 
of  predicting  certain  aspects  of  rotorcraft  response  is  still  poor. 

As  economic  development  continues  its  global  expansion,  the  world  market  is  be¬ 
coming  more  competitive.  Unfortunately,  both  the  production  and  operational  costs 
of  rotorcraft  have  been  too  high  and  have  made  the  rotorcraft  industry  noncompetitive 
in  the  commercial  transportation  marketplace,  Ref.  [1].  A  computer  based  simula¬ 
tion  tool  will  certainly  provide  a  solution  for  increased  productivity  and  reduced  cost. 
A  reliable  simulation  tool  will  provide  the  rotorcraft  industry  with  the  capability  to 
achieve  more  affordable  development,  test  and  operation  of  rotorcraft. 

1.1  Background 

Although  there  exists  a  number  of  tools  for  rotorcraft  modeling  and  analysis,  the 
prediction  of  the  helicopter  behavior  in  flight  including  in  the  design  phase  is  a  primary 
industrial  challenge.  Some  areas  that  require  model  enhancement  for  high  fidelity 
simulation  are  the  helicopter  off-axis  response  prediction,  aerodynamic  interference 
modeling,  and  simulation  of  unique  environments,  such  as  vortex  ring  state,  or  ship 
deck  landing. 
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1.1.1  Off- Axis  Response 

One  of  the  most  perplexing  problems  with  existing  rotorcraft  flight  analysis/simulation 
programs  is  the  inability  to  predict  the  correct  vehicle  off-axis  response  (e.g.,  lat¬ 
eral  axis  response  due  to  the  longitudinal  control).  Not  only  does  the  magnitude  not 
correlate  with  flight  test  data,  but  the  trend  is  also  erroneous,  Refs.  [2]  and  [3].  The 
current  inability  of  simulations  to  predict  the  correct  qualitative  rotorcraft  response 
limits  their  application  in  design,  flight  test  support  and  training. 

The  problem  has  attracted  worldwide  attention  and  efforts  in  resolving  it,  Refs.  [4] 
to  [8].  These  investigations  have  mainly  concentrated  on  the  pitch-roll  coupling.  They 
have  indicated  that  the  physics  which  causes  a  discrepancy  between  the  test  data  and 
the  simulation  in  the  off-axis  response  correlation  is  related  to  the  rotor  aerodynamics. 
Refs.  [4]  and  [5]  emphasize  the  dynamic  wake  distortion  during  transient  maneuvers 
of  the  rotorcraft.  When  the  rotor  tip-path-plane  tilts  due  to  the  application  of  the 
cyclic  control,  the  rotor  tip  vortex  deforms  in  a  dynamic  fashion  resulting  in  an 
induced  velocity  variation  over  the  rotor  disk.  The  influence  of  the  induced  flow  is 
such  that  it  causes  a  reversal  in  the  off-axis  tip-path-plane  tilt  (lateral  tilt  in  the  case 
of  pitch  and  longitudinal  tilt  in  the  case  of  roll)  in  a  direction  that  agrees  with  the 
flight  test  data. 

Reference  [4]  develops  a  model  based  on  a  prescribed  vortex  wake  trajectory  for 
steady  helicopter  body  rotation  in  hover.  Reference  [5]  modifies  the  inflow  dynamics 
from  simple  momentum  theory  in  hover  by  adding  a  constant  gain  feedback  term 
explicitly  related  to  the  helicopter  pitch  or  roll  rate.  A  separate  numerical  solution 
to  a  curved  vortex  tube  associated  with  the  steady  pitch  rate  is  applied  to  obtain 
the  parameters  for  the  inflow  dynamics  equations.  Reference  [6]  applies  parameter 
identification  techniques  to  flight  test  data  to  capture  this  effect  by  identifying  an 
aerodynamic  phase  angle  that  is  tuned  to  improve  the  off-axis  response  correlation 
with  the  test  data.  Reference  [7]  introduces  a  concept  of  steady  induced  angular 
momentum  due  to  the  rotor  downwash  swirl.  The  cross  coupling  moments  due  to 
the  so  called  “virtual  inertia  effect”  during  helicopter  pitch  and  roll  are  added  to  the 
helicopter  body  dynamics  equations.  The  moments  introduced  by  Ref.  [7]  behave  like 
mechanical  gyro  moments,  but  are  from  the  air  rotation.  Ref.  [7]  shows  that  adding 
these  moments  in  the  helicopter  body  dynamics  improves  the  pitch-roll  cross  coupling 
response.  However,  no  formulation  or  detailed  calculation  of  these  moments  has  been 
given  by  Ref.  [7].  Ref.  [8]  applies  the  vortex  method  using  nonuniform  vortex  rings 
in  a  dynamic  fashion  to  compute  the  induced  velocities  at  both  main  rotor  and  tail 
components.  Ref.  [8]  shows  a  better  prediction  in  phugoid  and  lateral-directional 
modes. 

The  above  efforts  have  greatly  enriched  our  understanding  of  the  cross  coupling 
response  problem.  The  investigations  have  indicated  that  the  pitch-roll  cross  coupling 
results  from  combined  effects  of  various  phenomena.  The  rotor  aerodynamics  is  the 
key  to  solving  the  poor  prediction  of  the  cross  coupling  response  for  existing  rotorcraft 
simulation  programs.  Although  some  simple  formulations  have  been  offered,  they  are 
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limited  to  steady  rotation  in  hover  and  rely  on  separate  methods  (system  identification 
or  vortex  method)  to  obtain  the  cross  coupling  parameters.  There  is  a  need  to  address 
this  problem  in  a  systematic  way  and  provide  a  well  defined  solution  that  is  valid  for 
the  full  range  of  the  flight  conditions  encountered  in  helicopter  flight  operation. 

1.1.2  Aerodynamic  Interference 

The  efforts  taken  in  Refs.  [4]  to  [7]  have  been  focused  on  main  rotor  pitch-roll  cross 
coupling.  There  are  a  broader  range  of  off-axis  response  phenomena  that  are  more 
closely  related  to  the  aerodynamic  interference.  Unfortunately,  this  is  another 
poorly  modeled  phenomena  in  rotorcraft  simulation. 

The  aerodynamic  interference  includes  interference  of  rotor-rotor,  rotor-fuselage, 
and  rotor-empennage  (vertical  fin  and  horizontal  surfaces).  These  aerodynamic  inter¬ 
ferences  can  have  considerable  influence  on  rotorcraft  performance,  flight  control  and 
vibration.  It  has  been  noticed  that  much  of  the  redesign  efforts  on  helicopters  over 
the  past  30  years  have  been  driven  by  the  unexpected  negative  impact  of  aerodynamic 
interference,  Refs.  [9]  and  [10].  The  aerodynamic  interaction  plays  a  major  role  in 
affecting  the  helicopter  off-axis  responses,  such  as,  pitch  due  to  sideslip,  yaw  due  to 
collective  or  pitch,  etc..  The  most  important  interference  effects  are  the  ones  from 
the  main  rotor  wake.  It  is  also  the  most  difficult  to  model  and  time  consuming  to 
compute.  There  have  been  both  analytical  and  experimental  works,  Refs.  [11]  to  [16] 
that  investigate  the  phenomena.  However,  at  the  present  time,  most  flight  simulation 
programs,  either  use  simple  momentum  estimation  models,  or  adopt  crude  empirical 
models  or  table  look-ups  for  specific  vehicles.  The  sophisticated  and  time  consuming 
vortex  wake  is  only  utilized  in  non-real  time  rotorcraft  trim  and  steady  state  flight 
analysis,  Refs.  [17]  to  [20]. 

1.1.3  Vortex  Ring  State 

Another  deficiency  of  existing  rotorcraft  simulations  is  flight  in  a  vortex  ring  state. 
Both  the  helicopter  main  rotor  and  tail  rotor  can  encounter  the  vortex  ring  state 
during  flight.  It  is  one  of  the  most  critical  flight  conditions  where  flight  safety  is 
concerned.  The  main  rotor  could  run  into  the  vortex  ring  state  when  the  helicopter  is 
descending.  The  tail  rotor  would  encounter  the  problem  when  it  is  in  sideward  flight 
or  in  a  side  wind.  The  capability  to  simulate  the  rotorcraft  in  the  vortex  ring  state 
will  provide  a  safety  guide  to  flight  test  and  training  that  will  reduce  the  hazards  of 
entering  the  vortex  ring  state  in  flight. 

The  vortex  ring  state  analysis  is  another  issue  that  needs  to  be  addressed  in  order 
for  a  full  flight  rotorcraft  simulation  program  to  be  applicable  in  all  ranges  of  the 
flight  environment.  Moreover,  a  reliable  vortex  ring  state  simulation  is  required  for 
flight  safety  in  support  of  flight  test  and  training.  There  are  experimental  studies 
describing  the  rotor  vortex  ring  phenomena,  Refs.  [21]  to  [26].  Empirical  uniform 
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induced  velocity  equations  for  vertical  descent  can  be  found  in  the  literature,  Ref. 
[27].  However,  very  limited  analytical  modeling  of  the  vortex  ring  state  is  available. 
Reference  [28]  presented  an  analytical  model  by  modifying  the  rotor  uniform  induced 
flow  in  the  vortex  ring  state  using  a  cylindrical  vortex  sheet  with  a  linear  decay  vortex 
strength.  The  model  is,  however,  limited  to  vertical  descent.  Recently,  ART  modified 
the  model  in  Ref.  [28],  allowing  the  equation  to  be  applicable  in  both  hover  and 
forward  flight,  Ref.  [16].  But,  the  modified  model  has  a  problem  in  transition  from 
normal  flight  to  the  vortex  ring  state. 

Another  interesting  aspect  of  the  vortex  ring  state  modeling  is  the  establishment  of 
the  vortex  ring  state  boundary.  The  establishment  of  the  vortex  ring  state  boundary 
can,  on  the  one  hand,  serve  as  a  safety  guide  for  rotorcraft  operation.  It  can  also 
be  used  to  modify  the  finite  state  dynamic  wake  equations  for  a  smooth  modeling 
transition  from  the  normal  flow  to  the  vortex  ring  state.  There  has  been  research  work 
performed  in  determining  the  vortex  ring  state  boundary.  Reference  [24]  assumed  that 
the  normal  component  of  the  blade  tip  vortex  translational  velocity  is  the  sum  of  half 
of  the  averaged  rotor  downwash  and  the  normal  relative  airflow  speed.  From  this 
assumption,  Ref.  [24]  then  defined  the  vortex  ring  state  boundary  by  setting  the 
normal  component  of  the  blade  tip  vortex  translational  velocity  to  zero.  Physically, 
this  boundary  criteria,  in  fact,  addresses  the  center  of  the  vortex  ring  state  instead  of 
the  boundary.  Reference  [29]  proposed  a  different  criteria  for  entrance  of  the  vortex 
ring  state  by  checking  if  the  component  of  the  relative  airflow  is  negative  in  the 
direction  of  the  rotor  downwash.  In  application,  this  criteria  is  too  conservative  since 
it  indicates  that  the  helicopter  will  enter  the  vortex  ring  state  during  the  descent 
at  almost  all  descending  angles.  In  reality,  the  helicopter  can  still  fly  safely  at  a 
small  descent  speed  without  entering  the  vortex  ring  state.  When  the  helicopter  is 
in  forward  flight,  the  safety  margin  is  even  larger.  Therefore,  a  more  accurate  vortex 
ring  state  criteria  is  yet  to  be  defined. 


1.1.4  Rotor  Wake  Modeling 

The  previously  identified  defects  in  current  rotorcraft  simulations  are  mostly  related 
to  the  rotor  aerodynamics.  Rotor  aerodynamics  remains  one  of  the  most  challeng¬ 
ing  tasks  in  rotorcraft  simulation  modeling.  Among  all  of  the  rotor  aerodynamic 
phenomena,  the  rotor  wake  problem  remains  the  essential  issue.  The  wake  is  of 
primary  importance  in  determining  the  rotor  aerodynamic  behavior,  which  directly 
affects  almost  every  aspect  of  rotorcraft  flight  including  performance,  flight  dynam¬ 
ics,  stability,  vibration,  and  aeroacoustic  noise.  Detailed  modeling  of  the  rotor  wake 
has  traditionally  relied  on  modeling  of  a  discretized  field  of  vortex  elements  that  are 
trailed  and  shed  from  the  rotor  blade  due  to  radial  and  azimuthal  changes  in  airloads. 
This  “vortex  wake”  can  predict  the  induced  velocity  distribution  from  the  potential 
functions  used  to  describe  each  discrete  vortex  element  in  the  model.  Two  types  of 
vortex  wake  models  are  commonly  used.  The  Prescribed  Wake  model  assumes  the 
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geometry  of  the  vortex  elements  relative  to  the  rotor  and  only  updates  the  strength 
of  these  elements  due  to  the  changing  airloads  on  the  rotor  blades.  The  Free  Wake 
model  allows  the  position  of  the  vortex  elements  to  change  due  to  the  local  induced 
velocity  at  each  element.  Both  approaches  are  iterative  and  enforce  a  periodic  steady 
state  condition  on  the  rotor  loads  to  compute  the  strength  and  position  of  the  vortex 
elements  at  each  flight  condition.  They  are  therefore  used  predominantly  for  steady 
state  performance  prediction. 

The  finite  state  dynamic  wake  model  developed  by  Peters  and  He,  Refs.  [30] 
to  [32],  represents  a  major  advance  in  rotorcraft  simulation  technology.  Beginning 
with  potential  functions  derived  from  a  prescribed  wake  geometry,  this  method  pos¬ 
tulates  a  truncated  series  solution  to  the  induced  velocity  distribution,  and  derives 
closed  form  equations  for  the  coefficients  of  this  solution.  The  level  of  detail  of  the 
flow  distribution  may  then  be  customized  to  the  user’s  application  by  the  selection 
of  the  number  of  radial  and  azimuthal  terms  in  the  expansion  (the  finite  state).  A 
key  feature  of  this  approach  is  the  inclusion  of  an  acceleration  potential  based  on  the 
inertia  of  the  air  mass  that  allows  for  the  prediction  of  the  transient  behavior  of  the 
inflow  (dynamic  inflow).  The  finite  state  dynamic  wake  model  has  the  sophistication 
of  the  classical  prescribed  vortex  wake  for  modeling  accuracy.  At  the  same  time  it  of¬ 
fers  efficiency  in  satisfying  the  stringent  computational  demands  of  flight  simulation. 
The  capability  of  the  finite  state  dynamic  inflow  model  to  predict  induced  velocity 
under  both  steady  and  maneuvering  flight  conditions  is  a  major  advantage  over  the 
vortex  wake  models  that  are  based  on  periodicity  assumptions  for  the  wake  and  are 
therefore  limited  to  steady  state  conditions.  The  increasingly  widespread  application 
of  finite  state  dynamic  inflow  approach  throughout  rotorcraft  flight  dynamics  com¬ 
munity  makes  it  an  ideal  candidate  for  enhancement  to  resolve  the  deficiencies  of 
existing  rotorcraft  simulation  models. 

1.1.5  Comprehensive  Interdisciplinary  Simulation 

A  distinctive  characteristic  of  rotorcraft  is  the  strong  coupling  of  its  structure,  aerody¬ 
namics,  propulsion,  and  control  system.  Any  meaningful  rotorcraft  simulation  must 
model  all  these  disciplines  and  their  couplings.  The  fidelity  of  the  simulation  highly 
depends  on  how  well  each  of  these  disciplines  and  their  couplings  are  addressed.  The 
FLIGHTLAB  simulation  program  developed  by  Advanced  Rotorcraft  Technology, 
Inc.  (ART)  offers  a  comprehensive,  multidisciplinary  rotorcraft  analysis/simulation 
software  tool,  Ref.  [33].  It  provides  physically  based  rotorcraft  modeling,  including 
models  of  different  levels  of  sophistication  (selective  fidelity)  for  the  structure,  aero¬ 
dynamics,  propulsion  (engine/drivetrain),  and  flight  control  system.  FLIGHTLAB 
models  are  reconfigurable  to  deal  with  different  types  of  rotorcraft  (single-main  rotor, 
tandem  rotors,  co-axial  rotors,  and  tiltrotors)  as  well  as  different  rotor  configurations 
(articulated,  hingeless,  bearingless,  teetering,  and  gimbaled).  The  reconfiguration 
can  be  done  by  selection  from  available  options  and  entering  different  data,  instead 
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of  modifying  the  simulation  code.  FLIGHTLAB  provides  a  unified  analysis  and  sim¬ 
ulation  environment  to  carry  out  the  performance,  stability,  control  response  (time 
and  frequency),  load  calculations,  and  pilot-in- the-loop  simulation.  The  enhancement 
of  rotorcraft  modeling  in  FLIGHTLAB  with  the  capability  of  predicting  off-axis  re¬ 
sponse,  modeling  aerodynamic  interaction  and  simulating  the  vortex  ring  state  will 
provide  a  valuable  computer  aided  simulation/ analysis  tool  to  support  flight  test  and 
flight  simulation  effectively. 


1.2  Objectives 

The  overall  objective  is  to  enhance  current  rotorcraft  simulation  modeling  to  offer 
a  cost  effective  tool  that  supports  land  and  sea  based  testing  and  operational  anal¬ 
ysis.  In  achieving  the  goal,  the  Phase  I  efforts  are  to  develop  a  theoretical  finite 
state  dynamic  wake  formulation  to  include  the  wake  distortion  effect  and  investigate 
the  feasibility  of  the  approach.  Under  Phase  I,  the  new  formulation  will  approach 
the  problem  in  a  systematic  way  so  that  the  resulting  model  will  be  applicable  in 
a  full  range  of  flight  conditions.  The  formulation  will  also  include  the  rotor  wake 
interference  for  aerodynamic  coupling  among  rotors  and  other  lifting  surfaces.  The 
formulation  will  be  developed  for  implementation  in  a  comprehensive  rotorcraft  sim¬ 
ulation  program  to  allow  for  the  interaction  of  various  disciplines  related  to  the  flight 
modeling  and  simulation,  such  as  the  coupling  of  the  finite  state  dynamic  wake  with 
the  unsteady  airloads,  the  structure  elasticity,  the  propulsion  system  dynamics  and 
flight  control  system. 


1.3  Approach 

The  approaches  taken  are  to  resolve  the  deficiency  of  current  rotorcraft  simulation 
programs  in  the  off-axis  response,  aerodynamic  interaction,  and  vortex  ring  state  pre¬ 
diction  through  the  improvement  of  the  Peters/He  finite  state  dynamic  wake  model, 
Refs.  [30]  to  [32].  The  enhanced  finite  state  dynamic  wake  model  will  model  the 
transient  wake  deformation  and  inherently  address  the  source  of  off-axis  correlation, 
while  providing  a  more  computationally  efficient  approach  than  a  vortex  wake  model. 
Also  the  enhancement  will  provide  a  more  global  and  robust  representation  than  a 
semi-empirical  model  tuned  to  the  experimental  data. 

The  enhanced  finite  state  dynamic  wake  model  will  be  incorporated  into  FLIGHT¬ 
LAB  to  evaluate  the  effectiveness  of  these  enhancements  in  improving  the  prediction 
of  the  flight  dynamic  response.  FLIGHTLAB’s  comprehensive  modeling  and  selective 
fidelity  capabilities  will  ensure  a  uniform  level  of  sophistication  in  modeling  al  inter¬ 
acting  phenomena.  For  Phase  One,  the  enhanced  formulation  for  the  general  finite 
state  dynamic  wake  model  will  be  developed.  The  implementation  is  made  to  a  three 
state  inflow  model  to  demonstrate  the  feasibility  of  the  approach.  A  full  implemen- 
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tation  of  the  general  finite  state  dynamic  wake  equations  will  be  addressed  in  Phase 
Two. 

The  details  of  the  formulation  of  the  rotorcraft  simulation  model  enhancement 
will  be  presented  in  the  chapters  that  follow.  The  innovations  of  this  approach  are 
first  summarized: 

1.  Inclusion  of  the  tip  vortex  distortion  effect  in  the  widely  used  finite  state  dy¬ 
namic  wake  model  to  overcome  the  poor  prediction  of  rotorcraft  off-axis  re¬ 
sponse. 

2.  Modification  of  finite  state  dynamic  wake  model  in  the  vortex  ring  state  to 
generate  a  simulation  model  applicable  across  the  flight  envelope. 

3.  Enhancement  for  the  finite  state  dynamic  wake  model  to  allow  for  aerodynamic 
interaction. 

4.  Integration  of  the  enhanced  finite  state  dynamic  wake  model  into  a  compre¬ 
hensive  simulation  program  (FLIGHTLAB)  to  address  the  inter-disciplinary 
couplings. 


Chapter  2 


Finite  State  Wake  with  Distorted 
Tip  Vortex 

Ref.  [34]  first  observed  that  the  theoretical  pitch-roll  damping  of  of  helicopters  did 
not  agree  with  flight  measurements  and  attributed  this  discrepancy  to  inflow  variation 
over  the  rotor  disk.  Ref.  [35]  provided  a  mathematical  model  to  explain  this  phe¬ 
nomenon,  and  made  the  first  systematic  exploration  that  established  a  quasi-steady 
relation  between  instantaneous  perturbation  in  thrust  and  perturbation  in  induced 
flow  based  on  momentum  theory.  Since  then,  the  theory  of  dynamic  inflow  has  expe¬ 
rienced  a  development  from  the  simple  momentum  approach  to  a  sophisticated  finite 
state  dynamic  wake  theory,  Refs.  [36]  to  [40]  and  Refs.  [30]  to  [32]. 

The  Peters/He  finite  state  dynamic  wake  model  describes  the  variation  of  the  rotor 
downwash  in  a  dynamic  fashion.  It  represents  the  induced  downwash  over  the  rotor 
plane  in  terms  of  a  series  of  polynomials  for  spanwise  variation  and  harmonic  (sine 
and  cosine)  functions  for  azimuthal  variation.  The  finite  state  dynamic  wake  model 
offers  a  formulation  so  that  the  user  can  tailor  the  sophistication  of  the  model  to  his 
needs  by  truncating  the  series.  It  is  applicable  in  both  steady  and  maneuver  flight 
conditions.  The  finite  state  dynamic  wake  model  presents  the  induced  flow  equations 
in  a  format  consistent  with  those  of  the  rotor  dynamics.  Moreover,  it  provides  closed 
form  solutions  for  its  coefficient  matrices,  which  significantly  reduces  the  computation 
time.  Because  of  these  advantages  over  other  wake  solutions,  the  finite  state  dynamic 
wake  model  has  found  widespread  applications  in  flight  dynamics  modeling  and  sim¬ 
ulation,  and  is  becoming  a  standard  industry  formulation  for  rotorcraft  simulation. 
The  original  Peters/He  finite  state  dynamic  wake  model,  Refs  [30]  to  [32],  has  three 
major  limitations.  First,  it  gives  the  analytical  closed  form  solution  to  the  normal 
downwash  at  the  rotor  plane  only,  which  limits  its  application  for  rotor  interference 
modeling.  Second,  its  induced  velocity  calculation  is  derived  from  a  prescribed  undis¬ 
torted  cylindrical  wake  geometry.  Third,  it  has  a  singularity  in  the  vortex  ring  state 
where  the  well  defined  slipstream  no  longer  exists.  This  prevents  it  from  accurately 
predicting  the  cross  coupling  response  and  the  vortex  ring  state. 
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2 . 1  Finite  State  Dynamic  Wake  Theory  of  Peters /He 

Before  we  proceed  to  the  derivation  of  an  enhanced  formulation  to  the  Peters/He  finite 
state  dynamic  wake  model  for  including  the  tip  vortex  distortion  effect,  a  brief  review 
of  the  fundamental  equations  of  the  finite  state  dynamic  wake  theory  is  necessary. 
The  Peters/He  finite  state  dynamic  wake  model  offers  a  new  approach  to  complicated 
helicopter  rotor  unsteady  wake  modeling.  It  solves  the  problem  from  fundamental 
unsteady  fluid  dynamics.  The  solution  is  based  on  the  governing  equations  for  an 
incompressible  potential  flow,  Ref.  [30]. 

For  an  incompressible  potential  flow  with  small  perturbations,  the  continuity  and 
momentum  equations  can  be  written  in  index  notation  as 

Qi,i  =  0  (2.1) 

Qi, o  -  VT  qi,s  =  i  (2.2) 

Eq.  (2.1)  is  the  continuity  equation  for  mass  conservation  and  Eq.  (2.2)  is  the 
momentum  equations  for  force  equilibrium.  The  above  equations  are  written  in  index 
notation  for  conciseness.  The  index  i  represents  i— component  of  the  velocity.  The 
are  the  perturbation  velocity  components,  <E>  is  the  pressure,  ( ),0  is  a  nondimensional 
time  derivative,  and  (  )^  is  the  derivative  along  the  relative  wind  direction. 

From  equation  (2.2),  it  can  be  seen  that  spatial  variation  of  the  pressure  is  the 
superposition  of  contributions  from  both  the  gradient  of  velocity  along  the  relative 
wind  direction  (or  convection)  and  the  unsteadiness  of  the  flow  field.  The  unsteady 
contribution  qhQ  is  unique  to  the  finite  state  induced  flow  theory.  All  conventional 
vortex  wake  models  do  not  include  this  term  and  therefore  are  only  quasi-steady  in 
nature. 

In  the  solution  of  Eqs  (  2.1  and  2.2),  we  divide  the  pressure  into  two  parts  (i.e., 
the  part  due  to  convection,  denoted  as  $y  ,  and  the  part  due  to  unsteadiness  denoted 
as  $x).  Then  the  pressure  can  be  expressed  as 


(2-3) 

or 

i,i  = 

(2.4) 

where 

o 

1 

II 

e 

(2.5) 

=  VTq^ 

(2.6) 

If  we  differentiate  equation  (2.2)  with  respect  to  the  index  i ,  then  interchange  the 
differentiation  order  of  both  (i,  0)  and  (£,  i )  with  the  help  of  continuity  equation  (2.1), 
a  Laplace’s  equation  for  pressure  function  can  be  obtained  as  follows, 

$  «  =  0 


(2.7) 
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In  a  similar  way,  from  equations  (2.5)  and  (2.6),  we  have 

=  0  (2.8) 

and 

=  0  (2-9) 

Equations  (2.8)  and  (2.9)  indicate  that  each  part  of  the  total  pressure  satisfies  Laplace’s 
equation.  Thus,  each  can  be  thought  of  as  an  acceleration  potential,  which  is  signifi¬ 
cant  for  the  formulation  of  the  theory.  The  boundary  conditions  for  each  $  compo¬ 
nent  are  that  the  pressure  match  the  blade  loading  on  the  rotor  blades  and  be  zero 

at  infinity. 

By  the  method  of  separation  of  variables,  the  solutions  to  the  Eq.  (2.7)  can  be 
obtained  as,  [30]. 

$  =  -i  f]  £  PnmMQ™(^)[(Cc)cosW)  +  (T™)sin(mi>)} 

^  771=0  71=771+1, 771+3, •" 


Similarly,  we  can  have  the  solution  for  the  Eqs.  (2.8)  and  (2.9). 

1  OO  0°  _  _ 

<tv  =  --Y.  E  +  (ov«nW<)] 

“  m=0  71=771+1, 771+3, ••• 

(2.11) 

-J  oo  o°  _ 

£  P^)Q^m(Tr)Acos(mi,)  +  K‘)Asin  (m>«] 

^  771=0  71=771+ 1,771+3,* •• 

(2.12) 

where  the  P™(  v)  and  Qn(ir})  are  the  normalized  first  and  second  kind  of  associated 
Legendre  functions.  The  t™c  and  t™3  are  pressure  functions  and  can  be  related  to 
the  blade  lift.  The  v  ,  rj  and  ^  are  the  ellipsoidal  coordinates  as  defined  in  Fig.  2.1. 

To  establish  the  inflow  dynamic  equations,  we  expand  the  normal  induced  flow  at 
rotor  plane  in  terms  of  radial  distribution  and  azimuthal  variation, 

OO  OO 

w(r,i>,t)  =  'jT  Z  C(r)[an(*)cosW0  +  Pn{t)sin(mip)}  (2.13) 

771=0  71=771+1, 771+3, ••• 


The  radial  expansion  functions,  0™(r),  are  the  polynomial  functions  in  the  radial 
position,  r.  The  ct™  and  are  the  induced  flow  expansion  coefficients.  They  are, 
in  fact,  the  states  of  the  rotor  induced  flow.  Once  these  inflow  states  are  solved  from 
the  inflow  dynamics  equations,  the  induced  flow  distribution  at  the  rotor  plane  can 
be  obtained  from  the  Eq.  (  2.13  ). 
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As  suggested  by  the  Eq.  (2.4),  the  induced  flow  field  can  be  solved  by  the  super¬ 
position  of  the  unsteady  and  the  convection.  Solving  the  Eq.  (2.5)  gives 


>]{<}'  =  \  {  rZ  J4 

(2.14) 

(2.15) 

At  the  same  time,  solving  the  Eq.  (2.6)  results 


{  <  ) 

_  1 
'  “  2 

;  V  ]  { if  }v 

(2.16) 

{  ' 

1  =  5 

[  V  ]  { rr  f 

(2.17) 

Combining  the  Eq.  (2.14)  with  (2.16)  and  the  Eq.  (2.15)  with  (  2.17)  and  using 
the  relationship  of  Eq.  (2.4)  gives  final  finite  state  induced  flow  equations  as 


[M]  {  <  }  +  [  ]"'  1  {  <  }  =  1 1  {  t?  } 


{  k  }  i  +[i* r i  i  w  i  [=;i  { } 


(2.18) 


(2.19) 


In  the  above,  the  r™c  and  r ™  are  the  inflow  forcing  functions.  They  can  be  computed 
from  blade  airloads,  Ref.  [32]. 


T“  =  if 

T"’  =  »CM*]«»(»>iM 


(2.20) 


(2.21) 


(2.22) 


where  Lq  is  blade  sectional  circulatory  lift  that  can  be  evaluated  from  any  time  domain 
unsteady  airload  model,  Ref.  [41]. 

In  the  above  equations,  the  [ L]c  and  [L]s  matrixes  are  obtained  from  a  nondis- 
torted  cylindrical  wake  geometry  in  closed  form.  It  does  account  for  the  effect  of  hub 
translational  motion  through  a  wake  skewing  angle,  Fig.  2.2.  But,  it  assumes  no 
wake  distortion  due  to  the  nonuniform  wake  induced  velocity  and  hub  rotation.  The 
enhancement  is  to  develop  a  formulation  that  includes  these  effects. 
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2.2  Formulation  with  Tip  Vortex  Distortion 

The  rotor  blade  tip  vortex  has  been  a  dominant  factor  that  affects  rotor  performance 
and  response.  Thus,  the  modeling  of  the  tip  vortex  has  long  been  a  central  issue  in 
the  rotor  solution.  The  existing  Peters/He  finite  state  wake  theory,  Refs.  [30]  to  [32], 
relies  on  adding  more  inflow  states  to  approximate  the  tip  vortex.  Since  the  tip  vortex 
is  a  concentrated  vortex  with  strong  discontinuity  while  the  inflow  states  are  continuos 
functions,  the  approximation  with  continuous  inflow  states  to  the  discontinuous  tip 
vortex  has  limitations.  The  approximation  could  be  significant  at  the  tip  region  if 
the  tip  vortex  is  strong,  such  as  in  the  case  of  a  rectangular  blade  tip. 

Since  the  blade  tip  vortex  is  a  dominant  phenomena  in  the  rotor  wake  modeling,  we 
focus  our  efforts  in  introducing  the  tip  vortex  distortion  into  the  new  finite  state  wake 
formulation.  In  order  to  achieve  computational  efficiency  while  retaining  necessary 
modeling  accuracy,  it  has  been  a  common  practice  in  most  comprehensive  rotorcraft 
vortex  wake  codes  that  only  the  tip  vortex  is  allowed  to  move  freely  for  wake  distortion 
modeling. 

In  this  approach,  we  will  derive  two  alternative  formulations.  The  first  general 
formulation  allows  for  a  combination  of  a  tip  vortex  with  the  finite  state  dynamic 
wake  equations  for  a  complete  dynamic  induced  flow  solution.  The  other  formula¬ 
tion  provides  parametric  equations  where  the  induced  flow  influence  coefficients  are 
extracted  utilizing  the  results  from  a  separate  tip  vortex  model. 

In  the  first  formulation,  the  existing  finite  state  dynamic  wake  model  is  modified 
and  combined  with  a  tip  vortex  model  in  one  module.  The  effects  of  tip  vortex  dis¬ 
tortion  both  global  (such  as  those  due  to  the  hub  rotation)  and  local  (such  as  those 
due  to  vortex  self-induced  distortion)  will  be  allowed  for  in  the  model.  This  approach 
has  the  advantage  of  a  more  accurate  description  of  the  physical  phenomena,  but 
will  be  very  computationally  intensive.  It  will  be  suitable  for  non-real  time  analysis 
on  more  complicated  problems,  such  as  blade-vortex  interaction,  vibratory  load  and 
aeroacoustic  noise.  The  parametric  formulation,  on  the  other  hand,  utilizes  predeter¬ 
mined  wake  distortion  influence  coefficients  and  adds  little  additional  computation 
the  original  finite  state  dynamic  wake  solution.  Thus,  the  parametric  formulation 
offers  computational  efficiency  and  satisfies  the  demand  for  real  time  simulation.  Due 
to  the  complexity  of  the  vortex  deformation,  only  the  global  geometric  wake  distor¬ 
tion  from  the  hub  rotation  and  quasi-steady  local  vortex  geometry  distortion  due  to 
the  self-induced  velocity  can  be  included  in  the  parametric  formulation. 

To  start  with,  we  modify  the  blade  tip  vorticity  modeling  by  introducing  a  con¬ 
centrated  tip  vortex.  The  nondistorted  cylindrical  wake  geometry  is  also  improved 
by  including  the  tip  vortex  distortion  due  to  both  the  hub  rotation  and  wake  induced 
velocity.  Since  the  current  finite  state  induced  flow  model,  Eqs  (2.18  to  2.22),  has 
already  accounted  for  the  contributions  from  the  blade  tip  vorticity  in  terms  of  con¬ 
tinuous  distribution  of  the  blade  tip  loading,  a  double  counting  of  these  effects  needs 
to  be  avoided. 
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Let  r™  be  the  induced  flow  forcing  function  computed  from  the  blade  lift  distri¬ 
bution,  Fig.  2.3.  We  define  a  new  inflow  forcing  function  {t™}up  that  has  only  the 
contribution  from  the  blade  tip  region  (dot  lines  filled  area),  Fig.  2.3.  By  modifying 
Eqs.  (2.14)  to  (2.17)  to  allow  the  induced  velocity  due  to  the  convection  to  be  com¬ 
puted  from  a  distorted  tip  vortex,  we  can  derive  a  new  set  of  induced  flow  equations 
that  include  tip  vortex  distortion  effects. 

For  unsteady  part,  Eqs.  (2.14  and  2.15),  we  have 
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For  the  convection  part,  Eqs.  (2.16  and  2.17), 
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where  the  {a™}tip  and  {(3™}tiP  are  tiie  contributions  from  the  tip  vortex.  They  can  be 
computed  from  the  induced  flow  distribution  over  the  rotor  plane  due  to  the  distorted 
tip  vortex.  Multiplying  both  sides  of  Eqs.  (2.25)  and  (2.26)  by  [Lc]-1  and  [Ls]  1, 
respectively,  results  in 
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The  final  finite  state  induced  flow  equations  with  the  tip  vortex  distortion  effect 
can  be  obtained  by  combining  Eq.  (2.27)  with  (2.23)  and  Eq.  (2.28)  with  (2.24). 


[m]{  r+ii']-1!  } 

=  |({  rr }  -  { rr  }  J  +  [  P  {  <  LP  +  5  {  c  It  (2-29> 

=  j({  T”‘  }  -  {  rr  }sp)  +  [  £•  ] ' "  {  K  }m  +  5  {  rr  It  (2.30) 

It  is  noted  that  the  above  equations  introduce  the  inflow  forcing  function  associ¬ 
ated  with  induced  flow  rate  change,  {r^c}^p  and  {r£w}$p.  In  steady  flight,  they  are 
negligible.  But  during  maneuver  flight,  we  need  to  compute  them  from  the  change 
of  induced  flow  associated  with  the  tip  loading.  To  do  this,  a  set  of  inflow  state 
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equations  for  the  tip  vortex  can  be  applied 

=5  <2-31> 

[  m  i  { w },; + m  “  i  *■  lp = 5  { i*  <2-32) 

From  the  solution  of  the  Eqs.  (2.31)  and  (2.32),  we  can  find 

{  =  2  [M}{<  };  (2.33) 

{  }t  =  2  [  M  1  (  &  )t  (2M) 

{TnC}tip  and  {t 'nC}tip  can  als°  be  obtained  directly  from  the  integration  of  the  blade 
tip  loading  that  is  associated  with  the  induced  flow  rate.  In  this  way,  no  additional 
inflow  state  equations  need  to  be  solved. 


2.3  Induced  Flow  due  to  Distorted  Tip  Vortex 

In  the  new  finite  state  induced  flow  formulation,  Eqs.  (2.29)  and  (  2.30),  the  new 
terms  {oc™c}tiP  and  {0™c}tiP  allow  for  the  effects  of  tip  vortex  distortion.  Let  the 
wtip(r ,  \j),  t)  be  the  induced  velocity  associated  with  the  distorted  tip  vortex,  from  Eq. 
(2.13)  we  have 

1  r2  7T  rl  _ 

K}«,  =  2^1  J0{TP»wUl,(r,ip,t)}drdi,  (2.35) 

{dn)tip  =  ~  /  [T'Pn(i/)cos(Tnip)wtiP(r,'ip,t)]drd'ip  (2.36) 

7 T  Jo  J 0 

1  /*27T  rl 

{Pn)tip  =  ~  /  /  {rP™(v)sin(rmp)wtip(r,  ip,  t)]drd-tp  (2.37) 

7T  Jo  JO 

Note  that  we  have  made  use  of  the  orthogonality  of  the  trigonometrical  and  associated 
Legendre  functions,  Ref.  [30],  in  deriving  the  above  equations. 

In  implementation,  we  can  compute  the  induced  velocity  directly  from  the  dis¬ 
torted  tip  vortex  using  a  vortex  model.  The  induced  flow  coefficients  {a™c}tiP  and 
{Pn°}tip  are  then  obtained  from  the  Eqs.  (2.35)  to  (  2.37)  and  added  to  the  final 
inflow  dynamics  equations,  Eqs.  (2.29)  and  (  2.30). 


2.4  Parametric  Formulation  of  the  Tip  Vortex  Ef¬ 
fects 

We  can  also  formulate  the  induced  flow  due  to  the  distorted  tip  vortex  in  parametric 
form  so  that  no  tip  vortex  has  to  be  explicitly  included  in  the  solution  of  finite  state 
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induced  flow  equations.  This  can  be  done  by  extracting  the  wake  distortion  influence 
coefficients  from  the  induced  velocity  resulting  from  a  separate  tip  vortex  model. 
Instead  of  using  the  total  contribution  from  the  distorted  tip  vortex,  it  is  much  more 
efficient  to  make  use  of  the  perturbation  induced  velocity  due  to  the  wake  distortion 
in  the  parametric  formulation.  To  do  so,  we  rearrange  Eqs.  (2.25  and  2.26), 


{ < } 

{«■} 


f  +  {  A<  } 

r+{  } 


(2.38) 

(2.39) 


where  {Aa™}  and  {A/3™}  are  the  perturbation  induced  velocity  due  to  the  wake 
distortion.  They  are 
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These  perturbation  inflow  due  to  the  wake  distortion  associated  with  the  blade  tip 
vortex  can  be  further  computed  as 
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where  [ALc]up  and  [AZd]tiP  are  the  wake  distortion  influence  coefficients  account¬ 
ing  for  the  perturbation  effects  of  the  local  tip  vortex  distortion  due  to  self-induced 
velocity.  [ALc]rot  and  [A Ls}TOt  are  the  wake  distortion  influence  coefficients  for  the 
perturbation  effects  of  the  global  geometric  distortion  of  the  tip  vortex  due  to  tip 
path  plane  rotation.  These  influence  coefficients  are  normalized  by  an  averaged  mo¬ 
mentum  theory  induced  flow  to  reduce  the  rotor  thrust  dependency  in  extracting  the 
coefficients. 

Substituting  Eqs.  (2.42)  and  (2.43)  into  (2.25)  and  (2.26)  and  combining  with 
Eqs.  (2.23)  and  (2.24),  we  obtain 
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where 


(2.44) 
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2.5  Tip  Vortex  Modeling 

Two  solution  methodologies  are  typically  adopted  in  fluid  dynamics:  Lagrangian  or 
Euler  methods.  The  Lagrangian  method  solves  the  time  variation  of  each  individual 
vortex  element,  while  the  Euler  method  deals  with  the  flow  field  point  instead  of 
following  the  individual  vortex  elements.  The  Lagrangian  approach  has  the  advantage 
of  being  able  to  trace  the  vortex  element  and  attain  a  time  accurate  wake  geometry 
solution.  This  is  particularly  important  for  maneuver  flight,  as  required  in  flight 
simulation  modeling.  The  Lagrangian  approach  applies  a  time  marching  scheme  that 
is  sometimes  susceptible  to  solution  instability  in  low  flight  speed.  The  instability 
can  be  both  physical  and  numerical.  As  observed  from  flow  visualization,  the  rotor 
wake  can  become  unstable  under  certain  low  speed  conditions.  On  the  other  hand, 
the  numerical  truncation  of  the  number  of  the  tip  vortex  elements  retained  in  the 
solution  can  also  cause  the  trajectory  of  the  vortex  elements  near  the  truncation 
boundary  to  diverge.  Extending  the  tip  vortex  far  down  the  stream  will  help  the 
convergence  of  the  solution.  For  computational  efficiency,  an  approximate  tip  vortex 
geometry  can  be  used  at  the  far  field  since  its  effect  on  the  rotor  is  less  significant 
as  compared  to  those  vortex  elements  that  are  in  the  vicinity  of  the  rotor  blades. 
Numerical  damping  is  also  often  applied  for  a  converged  solution  when  a  nonlinear 
vortex  solution  is  pursued,  but,  it  filters  out  wake  geometry  dynamics.  The  Euler 
approach  is  based  on  the  field  concept.  It  presents  the  flow  characteristics  in  term 
of  the  field  description.  For  rotorcraft  applications,  the  Euler  approach  allows  for 
enforcing  the  periodic  condition  under  steady-state  flight.  The  enforcement  of  the 
periodic  condition  through  numerical  relaxation  benefits  the  numerical  stability,  but 
its  application  is  limited  to  steady  flight  condition. 
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For  maneuver  flight  simulation,  the  rotor  periodic  condition  no  longer  exists  as 
the  pilot  constantly  applies  the  controls.  The  Lagrangian  approach  is  therefore  the 
only  viable  method  for  a  physically  realistic  solution. 

2.5.1  Tip  Vortex  Trajectory 

The  induced  velocity  field  in  the  vicinity  of  a  lifting  rotor  is  determined  by  the  rotor 
wake  geometry  and  wake  strength.  In  general,  it  is  a  nonlinear  phenomena  since  the 
wake  strength  and  geometry  interact  with  each  other.  On  one  hand,  the  change  in 
wake  strength  (e.g.,  circulation  of  the  vortex)  results  in  a  variation  in  the  associated 
induced  velocity.  On  the  other  hand,  the  induced  velocity  alters  the  wake  geometry. 
The  perturbation  in  the  wake  geometry  in  return  will  change  the  induced  velocity 
distribution  that  affects  the  rotor  airloads,  and  thus  the  wake  strength. 

In  general,  an  accurate  description  of  rotor  wake  geometry  can  only  be  obtained  by 
combining  the  effects  of  helicopter  rotor  hub  motion  (both  translational  and  rotation 
velocity),  rotor  self-induced  velocity,  wind,  and  interference  induced  velocity  from 
other  lifting  rotors  or  lifting  surfaces  (such  as,  wing,  fuselage,  and  etc.).  Due  to 
the  complexity  of  the  problem,  certain  approximations  have  to  be  made  in  order  to 
achieve  a  feasible  solution. 

In  the  conventional  free  vortex  model,  an  averaged  rotor  tip-path-plane  is  as¬ 
sumed  for  the  tip  vortex  trajectory  solution.  The  tip  vortex  trajectory  is  referred 
to  the  tip-path-plane.  The  relative  geometric  position  between  the  tip  vortex  and 
tip-path-plane  remains  unchanged  during  the  transition  from  one  rotor  revolution  to 
the  next.  A  typical  helicopter  transient  peak  response  is  less  than  about  two  to  three 
rotor  revolutions.  Therefore,  the  averaged  tip  vortex  geometry  is  no  longer  valid  for 
maneuver  flight  modeling.  To  overcome  this  hurdle,  we  take  an  approach  in  which 
the  tip  vortex  will  be  emanating  from  the  physical  blade  tip  at  each  time  instant  dur¬ 
ing  the  maneuver.  Thus,  the  tip  vortex  trajectory  follows  the  instantaneous 
blade  tip  response,  rather  than  the  averaged  tip-path-plane.  The  accurate 
tip  vortex  shedding  position  is  essential  to  catch  the  effect  of  hub  transient  transla¬ 
tion  and  rotation.  Once  the  tip  vortex  is  shed  into  the  flow  field,  it  will  move  with 
the  resultant  velocity  that  is  the  vector  superposition  of  the  wind,  rotor  self  induced 
velocity  and  interference  velocity  from  other  aerodynamic  lifting  surfaces. 

=  tU<  +  Vi  +  vM  (2.50) 

at 

where  x  is  the  vortex  element  moving  velocity,  vwind,  Vi  and  vint  are  the  wind, 
rotor  self-induced  and  interference  velocity  respectively. 

2.5.2  Prescribed  and  Free  Vortex 

There  are  two  approaches  that  are  commonly  utilized  in  rotor  wake  geometry  mod¬ 
eling.  One  is  the  prescribed  wake  method  in  which  the  rotor  wake  geometry  is  pre- 


18 


scribed.  The  wake  geometry  description  can  be  based  on  either  experimental  data 
(flow  visualization)  or  analytical  approximation,  such  as  using  averaged  induced  flow 
for  self-induced  effect  on  the  wake  geometry.  The  other  approach  is  the  free  vortex 
concept.  The  free  vortex  solution  allows  the  vortex  elements  to  move  in  the  flow  field 
resulting  from  the  superposition  of  wind  and  vortex  induced  velocity.  The  force-free 
condition  is  satisfied  when  the  vortex  segments  are  aligned  with  the  resultant  flow 
velocity.  The  free  vortex  solution  models  physical  distortion  of  the  rotor  wake.  But, 
there  are  hurdles  associated  with  the  rotor  free  vortex  approach.  Two  major  issues 
are  the  vortex  decay  and  the  vortex  interactions.  All  free  vortex  solutions  have  to 
rely  on  tuning  vortex  core  size  to  prevent  singularities  in  order  to  obtain  a  converged 
solution. 

2.5.3  Induced  Velocity  Calculation 

The  induced  velocity  is  computed  based  on  the  Biot-Savart  Law,  Ref.  [42].  The 
induced  velocity  at  a  flow  field  point  ‘p’  due  to  a  straight  vortex  line  with  constant 
circulation  is  calculated  as 

^  r  ..  _ 

vi  -  47rr°P  x  TbT> 

(rop  +  rapKl  ~ 

rbprap  ~  rdot  +  rc(fap  +  rbp~  2rdot) 

where  fap  and  rip  are  the  position  vectors  from  the  two  ends  of  the  vortex  line  (points 
‘a’  and  ‘b’  )  to  the  point  ‘p’.  The  scalar  rdot  is  defined  as  rdot  =  fap  •  f^,  rc  is  the  size 
of  the  vortex  core;  and  T  is  the  circulation  of  the  vortex. 

2.5.4  Global  Geometric  Distortion  Effect 

Most  existing  rotor  vortex  wake  models  do  not  include  the  effect  of  transient  rotor  hub 
motion  in  maneuver  flight.  The  finite  state  wake  model  does  let  the  rotor  wake  follow 
the  rotor  tip-path-plane  instantaneously  in  response  to  the  hub  translational  motion, 
but  the  rotor  wake  moves  as  a  whole  cylinder  without  the  time  step  effect  on  each 
individual  wake  element.  Moreover,  for  hub  rotation,  even  the  steady  motion  effect 
has  been  completely  ignored.  These  are  the  first  order  effects  that  must  be  addressed 
for  the  maneuvering  flight  simulation.  The  hub  translation  causes  the  rotor  wake  to 
be  skewed,  while  the  hub  rotation  produces  curved  wake  geometry  which  has  been 
found  to  be  significant  in  affecting  the  helicopter  pitch-roll  cross-axis  response.  In 
our  writing,  we  refer  to  the  hub  motion  induced  effects  as  global  geometric  wake 
distortion  effect.  The  rotor  tip  vortex  geometry  can  also  be  further  distorted  by  its 
self-induced  velocity  and  other  interference  velocities.  We  refer  to  that  as  the  local 
wake  distortion  effect.  The  local  distortion  effect  is  important  for  blade- vortex 
interaction,  aeroacoustics,  and  vibrations. 


(2.51) 

(2.52) 
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The  first  step  we  take  for  enhancement  of  rotorcraft  simulation  modeling  is  to 
include  the  global  geometric  wake  distortion  effect  in  the  finite  state  wake  equations. 
To  model  the  effect  properly,  an  accurate  description  of  the  rotor  hub  motion  is 
essential.  The  comprehensive  rotorcraft  simulation  program — FLIGHTLAB,  Refs. 
[33],  [41],  [43],  and  [44].  provides  an  ideal  tool  due  to  its  geometrically  exact  nonlinear 
structural  dynamics  modeling.  All  the  airframe  dynamic  modeling  in  FLIGHTLAB 
is  fully  nonlinear  without  ignoring  any  higher  order  geometric  terms.  This  is  of 
importance  under  maneuver  conditions  where  the  attitude  change  of  the  helicopter 
body  can  be  large  enough  so  that  the  small  angle  assumption  is  no  longer  valid. 

In  our  tip  vortex  geometry  calculation,  the  tip  vortex  trajectory  is  described  rela¬ 
tive  to  the  physical  rotor  hub  frames  that  experience  linear  and  angular  motion  in  all 
three  directions.  The  rotor  hub  frame  is  attached  at  the  center  of  the  rotor  hub.  It  is 
defined  such  that  ‘x’  points  downstream,  ‘y’  t0  the  right  and  ‘z’  is  up.  The  relative 
position  of  the  tip  vortex  to  the  rotor  hub  frame  is  computed  at  each  time  step,  and 
thus  includes  the  instantaneous  change  in  the  rotor  hub  motion  under  maneuver  flight 
conditions.  In  the  model,  the  tip  vortex  emanates  from  the  physical  blade  tip  and 
then  sets  off  to  move  in  the  flow  field  that  results  from  the  wind  and  induced  veloc¬ 
ity.  To  extract  the  global  geometric  effect  on  the  wake  distortion,  the  vortex  moving 
velocity  is  assumed  to  be  that  of  rotor  momentum  theory  value.  Under  unsteady 
flight  conditions,  the  momentum  theory  value  is  updated  at  each  time  step  from  the 
instantaneous  rotor  thrust. 


2.5.5  Local  Distortion  Effect 

The  local  distortion  effect  here  refers  to  the  wake  distortion  due  to  the  nonuniform 
induced  velocity  field.  Including  this  effect  produces  the  so  called  free  vortex  wake 
model.  The  local  wake  distortion  effect  has  been  found  to  have  an  effect  on  the  rotor 
vibratory  loads  and  aeroacoustics.  Its  impact  on  the  vehicle  flight  dynamics  remains 
to  be  investigated.  In  Phase  One,  we  concentrate  on  the  finite  state  wake  modeling 
with  the  global  geometric  wake  distortion  effect.  We  will  proceed  to  model  the  local 
distortion  phenomena  and  investigate  its  effect  on  the  flight  simulation  under  Phase 
Two  effort. 

2.5.6  Vortex  Core  Model 

The  vortex  core  model  is  always  a  challenging  issue  in  a  distorted  vortex  solution. 
The  core  size  of  the  vortex  plays  an  important  role  when  the  vortex  element  is  close  to 
the  point  where  the  induced  velocity  is  to  be  calculated.  Too  small  a  vortex  core  size 
can  cause  an  unrealistically  large  fluctuation  of  the  induced  flow  velocity  and  often 
lead  to  diverged  solutions  when  it  couples  with  rotor  dynamics.  On  the  other  hand,  if 
the  core  size  is  too  large,  an  inaccurate  solution  may  result.  Many  existing  free  wake 
codes  adopt  an  empirical  number  and  fix  the  core  size  for  the  whole  solution  process. 
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Instead  of  using  a  fixed  core  size  during  the  process  of  a  free  vortex  solution,  in 
this  approach,  we  allow  the  size  of  the  vortex  core  to  vary  with  the  time.  The  formula 
is  based  on  the  classical  incompressible  viscous  flow  solution.  In  this  solution,  a  two 
dimensional  plate  is  set  to  move  at  a  constant  forward  speed,  the  vortex  created 
from  the  plate  boundary  layer  is  then  grown  through  the  flow  field.  The  vortex  core 
variation  with  time  is  thus  obtained  as 


rvc  =  T]cVvet  (t  <ti) 

=  rvc o  +  rjcVuet  (t  >  tx)  (2.53) 

where  rjc  is  a  constant,  v  is  the  kinematic  viscosity  and  e  is  the  eddy  viscosity  coef¬ 
ficient.  The  parameter  tx  is  the  time  when  the  vortex  passes  by  the  first  following 
blade.  The  parameter  rvc0  is  the  core  size  change  at  the  time  of  tx.  This  formu¬ 
lation  allows  for  the  effect  of  blade/vortex  interaction  on  the  growth  of  the  vortex 
core.  When  applying  the  formula  to  the  helicopter  blade  tip  vortex,  the  physical  tip 
vortex  development  is  more  complicated  than  the  decay  process  of  the  ideal  laminar 
vortex  where  the  eddy  viscosity  coefficient  is  unity.  Thus,  the  parameter  e  needs  to 
be  selected  based  on  experimental  data. 

2.6  Results  and  Discussion 

The  results  presented  in  this  section  include  both  the  induced  downwash  variation 
due  to  the  tip  vortex  distortion  and  its  effects  on  the  helicopter  off-axis  response.  The 
inflow  influence  coefficients  extracted  from  the  distorted  tip  vortex  are  also  discussed. 


2.6.1  Wake  Distortion  Influence  Coefficients 

Figure  2.4  shows  the  trajectory  of  the  rotor  tip  vortex  under  three  different  tip-path- 
plane  pitch  rates  (  qtpp  —  0, 5.0, 10.0, 15.0  deg/ sec  )  during  hub  rotation.  Eight  rotor 
wake  revolutions  are  used  in  these  models.  The  wake  geometries  are  plotted  8  rotor 
revolutions  after  the  start  of  the  hub  rotation.  At  this  time,  the  hub  rotation  effect  has 
been  transported  to  the  last  wake  element.  The  tip  vortex  trajectories  that  emanate 
from  the  physical  blade  tip  present  two  distinctive  features.  The  rotation  of  the  rotor 
hub  produces  a  curved  vortex  geometry  as  compared  to  a  straight  cylindrical  line 
under  the  conventional  method.  Moreover,  the  hub  rotation  causes  the  vortex  sheet 
to  compress  on  one  side  and  expand  on  the  other.  The  extent  to  which  the  tip  vortex 
geometry  is  curved  and  compressed/expanded  depends  on  the  rate  at  which  the  hub 
is  rotating. 

Figure  2.5  shows  the  wake  geometries  at  30  knots  with  the  four  tip-path-plane 
rotation  rate.  The  curvature  of  the  wake  is  less  pronounced  than  that  in  hover, 
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however,  the  wake  concentration  near  the  rotor  trailing  edge  is  more  than  in  the 
similar  region  shown  in  Fig.  2.4. 

In  the  parametric  formulation,  the  elements  of  [ALc]rot  and  [ALs}TOt  are  the  global 
geometric  wake  distortion  influence  coefficients  due  to  the  hub  rotation.  In  modeling, 
the  blade  tip  vortex  is  trailed  from  physical  blade  tip  during  the  hub  rotation.  Thus, 
the  wake  distortion  effect  is,  in  fact,  directly  related  to  the  tip-path-plane  rotation 
as  formulated  in  Eqs.  2.42  and  2.43.  To  extract  the  influence  coefficients  from  the 
distorted  vortex  results,  it  needs  to  identify  the  independent  variables  with  which 
these  influence  coefficients  vary. 

Fig.  2.6  shows  the  variation  of  wake  distortion  influence  coefficients  during  hub 
pitch  and  roll  in  hover.  The  AT12  is  the  fore-to-aft  inflow  state  )  change  due  to 
the  pitch  rate.  The  AL2i  is  the  side-to-side  inflow  state  )  change  due  to  the  roll 
rate.  The  A Lu  and  A L22  are  the  cross  coupling  for  the  change  of  a]  due  to  the  roll 
and  the  due  to  the  pitch,  respectively.  It  is  noticed  that  the  influence  coefficients 
are  essentially  constants  soon  after  reaching  the  limiting  value,  i.e.,  independent  of 
pitch  angle.  It  is  also  noted  that  the  cross  coupling  between  the  induced  flow  change 
and  pitch /roll  motion,  i.e.,  the  A Ln  and  A L22  are  muchTess  as  compared  to  the 
on-axis  inflow  change  (AZq2  and  AI2i).  The  wake  distortion  influence  coefficients 
obtained  here  are  based  on  the  formulation  of  Eqs.  2.44  and  2.45.  If  we  convert  the 
influence  coefficient  to  an  equivalent  parameter  in  Ref.  [5],  we  have  an  equivalent 
value  of  1.45  as  compared  to  the  1.5  reported  in  the  Ref.  [5].  This  is  due  to  the 
fact  that  we  have  used  physically  more  accurate  helix  vortex  segments  instead  of  the 
vortex  ring  tube  as  adopted  in  Ref.  [5]. 

Fig.  2.7  shows  the  perturbation  induced  flow  state  variation  with  tip-path-plane 
angle  of  attack  (i.e.,  the  angle  between  the  free  stream  and  tip-path-plane)  during  the 
pitching  motion  at  different  pitch  rates  in  10  Knots  forward  flight.  The  inflow  states 
have  been  normalized  by  the  averaged  momentum  uniform  inflow  value.  Fig.  2.8 
shows  similar  results  at  80  Knots  forward  flight.  One  distinguished  feature  of  these 
results  is  the  linear  variation  with  the  pitch  rate  as  indicated  by  a  nearly  constant 
factor  difference  in  the  results  for  different  pitch  rates.  This  is  further  confirmed 
by  directly  plotting  of  the  influence  coefficients,  i.e.,  the  slopes  of  the  inflow  states 
with  respect  to  the  rotation  rate,  against  the  tip-path-plane  angle  of  attack,  Figs. 
2.9  and  2.10.  Fig.  2.9  shows  the  influence  coefficients  of  fore-to-aft  induced  flow 
distribution  due  to  pitch  rate  at  different  forward  flight  speeds.  The  variations  of 
influence  coefficients  with  various  pitch  rates  fall  into  essentially  one  curve  at  each  of 
the  forward  flight  speed  (10,  30,  60  and  80  Knots).  This  suggests  that  these  influence 
coefficients  can  be  either  tabulated  or  curve  fit  as  a  function  of  the  flight  speed  and 
the  tip-path-plane  angle  of  attack.  Fig.  2.10  presents  the  influence  coefficients  of 
side-to-side  induced  flow  distribution  due  to  pitch  rate.  This  is  the  inflow  due  to 
cross-axis  rotation.  As  seen,  the  effect  of  the  cross  coupling  is  one  order  smaller  as 
compared  with  the  on-axis  results  in  Fig.  2.9. 

Fig.  2.11  shows  the  influence  coefficients  of  side-to-side  induced  flow  distribution 
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due  to  the  roll  at  various  flight  speeds.  These  are  the  inflow  influence  coefficients  from 
the  on-axis  rotation.  As  shown,  the  influence  coefficients  reach  a  nearly  constant  value 
for  all  the  speeds,  indicating  that  the  side-to-side  induced  flow  state  is  only  a  linear 
function  of  the  roll  rate.  Fig.  2.12  presents  the  influence  coefficients  of  fore-to-aft 
induced  flow  distribution  due  to  the  roll  at  various  flight  speeds.  These  are  the  inflow 
influence  coefficients  from  the  cross-axis  coupling.  They  are  one  order  smaller  than 
the  influence  coefficients  for  the  on-axis  rotation.  Similarly,  the  influence  coefficients 
reach  nearly  constant  limiting  values  at  each  of  the  flight  speeds  and  are  basically 
independent  of  the  roll  rate. 

For  implementation,  the  inflow  influence  coefficients  due  to  the  wake  distortion 
are  also  plotted  against  forward  flight  speed,  Figs.  2.13  to  2.16.  Fig.  2.13  is  the 
variation  of  fore-to-aft  inflow  influence  coefficients  due  to  pitch  with  the  flight  speed. 
The  coefficients  are  plotted  for  three  typical  tip-path-plane  angle  of  attack.  A  large 
variation  occurs  in  the  speed  range  of  hover  to  30  Knots.  The  wake  distortion  effect 
decreases  with  the  flight  speed.  It  becomes  negligible  after  80  Knots.  Fig.  2.14  shows 
the  variation  of  side-to-side  inflow  influence  coefficients  due  to  roll  with  the  flight 
speed.  The  effect  of  the  wake  distortion  on  the  side-to-side  inflow  influence  coefficient 
decreases  as  helicopter  flies  forward.  Similarly,  the  effect  is  negligible  after  80  Knots. 
Figs.  2.15  and  2.16  are  the  variations  of  cross-axis  inflow  influence  coefficients  with 
forward  flight  speed.  The  results  show  that  these  cross-axis  inflow  influence  coeffi¬ 
cients  are  very  small  at  the  typical  range  of  tip-path-plane  rotation  (5,  10,  and  20 
degrees).  Thus,  the  effects  of  the  wake  distortion  on  the  cross-axis  inflow  influence 
coefficients  is  essentially  negligible. 

The  above  results  present  a  breakthrough  in  modeling  the  global  geometric  wake 
distortion  effect  in  the  finite  state  dynamic  wake  formulation  in  forward  flight  because 
all  the  work  before  this  have  been  limited  to  hover  flight.  Thus  far,  no  systematic 
results  have  been  obtained  for  the  wake  distortion  influence  coefficients  in  forward 
flight.  With  the  above  results,  we  can  apply  these  wake  distortion  influence  coefficients 
in  the  parametric  formulation  for  flight  simulation.  Thus,  our  accomplishments  have 
presented  a  practically  feasible  solution  to  the  problem. 

To  summarize,  the  perturbation  induced  flow  from  the  wake  distortion  due  to  the 
hub  rotation  is  a  function  of  flight  speed  and  tip-path-plane  angle  of  attack  for  pitch, 
and  flight  speed  and  tip-path-plane  lateral  tilt  angle  for  roll.  Most  importantly,  they 
are  only  linear  functions  of  the  pitch  and  the  roll  rates.  Thus,  there  is  no  hurdle  in 
implementing  the  parametric  formulation  to  allow  for  global  geometric  wake  distortion 
in  flight  simulation. 

2.6.2  Inflow  State  Variation  with  Vortex  Wake  Local  Distor¬ 
tion 

In  this  section,  we  examine  the  effect  of  vortex  local  distortion  due  to  rotor  self- 
induced  velocity  by  comparing  the  inflow  state  variation  between  the  free  vortex  and 
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the  finite  state  dynamic  wake  models.  The  inflow  states  for  the  vortex  wake  model 
are  extracted  from  the  induced  flow  distribution  over  the  rotor  plane  by  using  Eqs. 
2.35  to  2.37. 

Fig.  2.17  illustrates  the  free  vortex  wake  geometry  of  a  UH-60  rotor  in  30  Knots 
forward  flight.  This  is  the  speed  range  where  the  rotor  wake  distortion  is  the  most 
significant.  Fig.  2.18  compares  the  uniform  and  first  harmonic  inflow  states  com¬ 
puted  from  undistorted  and  distorted  wake.  These  are  the  three  most  important 
inflow  states  concerning  the  flight  dynamics.  They  are  the  inflow  states  describing 
uniform  (a°  ),  fore-to-aft  (a\  )  and  side-to-side  )  distributions.  As  shown,  both 
the  free  vortex  model  and  finite  state  dynamic  wake  results  have  shown  the  same 
variation  trend.  There  are  only  magnitude  difference  that  could  be  compensated  by 
the  trim  setting.  Thus,  the  effect  from  the  local  wake  distortion  will  be  higher  order 
as  compared  to  the  global  geometric  wake  distortion  that  affects  even  the  direction 
in  which  the  rotorcraft  responds  in  maneuver  flight.  Fig.  2.19  compares  the  inflow 
states  of  (a\  )  and  )  in  forward  flight.  The  oscillation  of  the  inflow  states  is  mainly 
due  to  the  blade  passage  effect  (or  blade  bound  vortex  effect).  The  4/rev  oscillation 
of  the  inflow  states  for  the  4-bladed  rotor  is  clearly  seen  for  both  the  undistorted 
and  distorted  wake.  The  wake  distortion  effect  is  not  significant  on  the  oscillation 
magnitude.  There  is  a  minor  phase  shift  for  the  inflow  states  that  may  be  due  to  the 
difference  in  inflow  collocation  points.  When  we  extract  the  wake  distortion  influence 
coefficients,  these  bound  vortex  contribution  will  not  be  a  concern  since  there  is  no 
bound  vortex  geometry  distortion. 

2.6.3  Off-Axis  Response 

The  results  presented  in  this  section  are  to  validate  the  enhanced  finite  state  dy¬ 
namic  wake  model  as  described  above.  The  UH-60  helicopter  response  is  computed 
in  FLIGHTLAB  simulation  environment.  The  rotor  blades  are  modeled  as  elastic 
using  the  FLIGHTLAB  modal  component.  The  enhanced  finite  state  dynamic  wake 
equation  is  applied  for  the  rotor  inflow  dynamics  modeling. 

The  first  case  investigated  is  the  UH-60  response  to  a  step  lateral  stick  in  hover. 
Figure  2.20  displays  the  pilot  control  inputs  as  applied  to  the  simulation.  Figure  2.21 
is  the  helicopter  body  response  to  the  pilot  stick  inputs.  The  results  from  the  model 
with  and  without  the  wake  distortion  due  to  the  hub  rotation  are  compared  to  the 
flight  test  measured  data.  As  seen,  the  model  without  considering  the  wake  geometry 
change  due  to  the  hub  rotation  predicts  a  pitch  rate  (off-axis)  response  opposite  to  the 
test  data.  The  wrong  prediction  of  the  response  trend  is,  however,  corrected  by  the 
enhanced  3-state  inflow  model  that  includes  the  wake  geometry  distortion  effect  due 
to  the  hub  rotation  during  the  transient  maneuver.  Figure  2.22  is  the  response  of  the 
helicopter  attitude.  The  enhanced  model  predicts  the  right  pitch  attitude  variation 
due  to  the  lateral  stick  input. 

To  understand  the  improved  prediction  from  the  enhanced  model,  Fig.  2.23  com- 


24 


pares  the  rotor  side-to-side  induced  flow  distribution  and  rotor  hub  pitch  moment 
from  the  two  models.  The  upper  plot  displays  the  induced  flow  distribution  at  the 
peak  of  the  body  roll  rate  after  the  pilot  lateral  stick  is  applied.  Compared  to  the 
existing  finite  state  inflow  model  without  the  hub  rotation  effect,  the  enhanced  model 
that  introduces  the  hub  rotation  effect  on  the  inflow  calculation  predicts  a  dramatic 
change  in  the  rotor  side-to-side  induced  flow  variation.  The  change  in  the  side-to-side 
inflow  variation  directly  affects  the  rotor  load  distribution  that  generates  a  much  dif¬ 
ferent  rotor  pitch  moment  response,  as  indicated  in  the  lower  plot  of  Fig.  2.23.  The 
hub  pitch  moment  change  from  the  enhanced  model  results  in  a  correct  helicopter 
pitch  response  to  the  lateral  stick  input. 

Using  the  enhanced  formulation,  we  have  further  investigated  the  effects  on  off- 
axis  response  in  forward  flight.  Figure  2.24  shows  the  UH-60  helicopter  response  to 
a  one  inch  aft  longitudinal  step  stick  input  at  60  Knots  forward  flight.  The  results 
from  three  models  and  their  correlation  with  the  flight  test  data  are  compared.  The 
first  simulation  model  uses  the  unsteady  uniform  inflow  with  Glauert  fore-to-aft  inflow 
distribution.  The  other  two  are  the  3-state  inflow  models  with  and  without  tip  vortex 
distortion.  For  the  on-axis  pitch  response  correlation,  all  the  three  models  have 
predicted  the  right  trend.  The  uniform  inflow  model  slightly  overshoots.  The  3-state 
inflow  model  without  the  wake  distortion  slightly  underestimates  the  response  peak. 
The  wake  distortion  model  provides  the  best  correlation.  But,  all  the  three  models 
matched  the  measured  data  within  a  reasonable  accuracy.  The  off-axis  response 
comparisons,  however,  have  shown  large  difference.  The  uniform  inflow  model  fails  in 
predicting  even  the  right  trend.  The  3-state  inflow  without  tip  vortex  distortion  shows 
too  small  of  a  roll  rate  amplitude  response.  Inclusion  of  the  tip  vortex  distortion  in 
the  3-state  inflow  model  not  only  predicts  the  right  direction  in  which  the  helicopter 
rolls,  but  also  improves  the  correlation  of  the  amplitude  of  the  roll  rate  response, 
although  there  is  a  remaining  difference  as  compared  to  the  flight  test  data.  Fig.  2.25 
shows  the  rotor  tip-path-plane  variation  during  the  maneuvering.  A  larger  lateral 
tilt  of  the  tip-path-plane  to  the  right  occurs  with  the  simulation  that  models  the 
wake  distortion.  The  wake  distortion  effect  is  also  reflected  on  the  rotor  hub  moment 
variation,  Fig.  2.26. 

It  is  noted  that  the  response  peak  of  the  roll  rate  predicted  is  still  less  than  that 
of  measured  data.  In  fact,  this  is  no  surprise  since  there  are  other  factors  that  affect 
the  helicopter  response  in  forward  flight.  This  remaining  difference  can  be  due  to 
the  fuselage  aerodynamic  interaction  on  the  rotor.  As  shown  in  Fig.  2.27  (from  Ref. 
[30])  there  is  a  rotor  flow  field  variation  due  to  the  existence  of  a  fuselage.  In  forward 
flight,  the  fuselage  induces  upwash  in  the  forward  part  of  the  rotor  plane  and  increases 
the  downwash  at  the  rear  portion  of  the  rotor  plane.  The  rotor  will  respond  to  this 
fore-to-aft  inflow  variation  in  the  lateral  flapping  that  will  further  increase  the  roll 
response  of  the  helicopter.  The  fuselage  beneath  the  rotor  mainly  affects  this  fore-to- 
aft  rotor  inflow  variation.  It  also  increases  with  forward  flight  speed.  In  Phase  Two 
effort,  we  will  develop  the  fuselage  interference  model  and  investigate  the  fuselage 


upwash  influence  on  the  helicopter  response. 


26 


*w 


Figure  2.1:  Ellipsoidal  Coordinates 
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Figure  2.6:  Global  geometric  wake  distortion  influence  coefficient 
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Figure  2.7:  Variation  of  normalized  induced  flow  states  with  Tpp  angle  of  attack  at 
10  Knots 
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Figure  2.8:  Variation  of  normalized  induced  flow  states  with  Tpp  angle  of  attack  at 
80  Knots 
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Figure  2.9:  Fore-to-aft  inflow  influence  coefficient  due  to  pitch  in  forward  flight 
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Figure  2.10:  Side-to-side  inflow  influence  coefficient  due  to  pitch  in  forward  flight 
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Figure  2.11:  Side-to-side  inflow  influence  coefficient  due  to  roll  in  forward  flight 
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Figure  2.17: 
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Figure  2.19:  Oscillating  inflow  states  from  distorted  and  undistorted  wake,  30  Knots 
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Figure  2.20:  Pilot  control  stick  inputs 
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Figure  2.24:  Wake  distortion  effect  on  off-axis  response  in  60  Knots 
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Chapter  3 

Finite  State  Inflow  Interference 


Including  the  tip  vortex  distortion  in  finite  state  dynamic  wake  formulation  improves 
the  prediction  of  pitch-roll  cross  axis  coupling.  There  is  another  class  of  off-axis  cou¬ 
pling  phenomena  associated  with  aerodynamic  interference.  For  example,  there  is  a 
pitch  response  from  sideslip  due  to  the  asymmetric  distribution  of  main  rotor  down- 
wash  over  the  horizontal  tail  surfaces.  Unfortunately,  the  aerodynamic  interference  is 
another  poorly  modeled  phenomena  in  rotorcraft  simulations.  Most  flight  simulation 
programs  at  the  present  time  use  simple  momentum  estimation  models,  or  rely  on 
crude  empirical  models  or  table  look-ups  for  specific  vehicles. 

The  simple  models  do  not  provide  enough  fidelity  to  generate  an  acceptable  sim¬ 
ulation.  The  empirical  models  or  table  look-up  methods  are  limited  to  a  specific  he¬ 
licopter  and  can’t  include  transient  dynamic  effects.  The  vortex  wake  models,  Refs. 
[17]  to  [20],  on  the  other  hand,  are  time  consuming  and  have  limitations  in  modeling 
maneuvering  flight.  Calculation  of  the  aerodynamic  interference  from  the  finite  state 
dynamic  wake  model,  Refs.  [30]  to  [32],  provides  a  solution  to  the  problems.  The  fi¬ 
nite  state  dynamic  wake  model  has  the  sophistication  of  the  vortex  wake  for  modeling 
accuracy  and  is  applicable  in  transient  flight.  At  the  same  time,  it  offers  efficiency 
in  satisfying  the  stringent  computational  demand  of  flight  simulation.  Moreover,  the 
finite-state  modeling  of  complicated  rotor  wake  dynamics  presents  a  solution  format 
that  is  consistent  with  all  existing  flight  dynamics  formulations.  The  finite-state  in¬ 
duced  flow  model  was  first  developed  to  compute  normal  induced  velocity  at  the  rotor 
plane.  But,  the  inflow  solution  has  been  limited  at  the  rotor  plane  and  can’t  be  used 
for  interference  modeling.  A  natural  and  beneficial  step  is  to  enhance  the  rotor  finite 
state  wake  modeling  capability  to  calculate  the  rotor  induced  velocity  at  any  arbi¬ 
trary  point  in  the  flowfield.  This  provides  solutions  for  main  rotor  interference  effects 
on  the  tail  rotor,  fuselage,  and  tail  aerodynamic  surfaces  (vertical  fin  and  horizontal 
surface).  In  a  previous  effort,  we  have  developed  a  basic  formulation  for  the  finite 
state  rotor  interference  solution.  Here  the  rotor  inflow  interference  model  is  further 
made  to  work  under  the  finite  state  dynamic  wake  model  enhanced  with  the  wake 
distortion.  We  have  also  performed  validation  of  the  finite  state  rotor  interference 
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model  and  investigated  its  effects  on  rotorcraft  trim  and  response. 


3.1  Formulation 


From  the  basic  equations  presented  in  the  previous  section,  the  induced  velocity  at 
an  arbitrary  flow  field  point  can  be  attained  through  an  integration  of  Eq.  (2.6).  This 
gives 

{3A) 

This  relationship  is  derived  in  the  relative-wind  frame  (a^,  y?,  z?),  Figure  3.1, 
but  it  is  also  true  in  the  wind-hub  reference  frame  (xw,  yW)  zw)  because  of  the  linear 
transformation  relationship  between  the  two  frames, 
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3.1.1  Partial  Derivative  Relationship 

Since  the  solution  of  the  pressure  function  $v  is  expressed  in  ellipsoidal  coordinates, 
Fig.  2.1,  the  calculation  of  induced  velocity  by  Eq.  (3.1)  needs  an  application  of 
the  chain  rule  for  the  evaluation  of  In  the  following,  a  matrix  of  first  partial 
derivatives  between  the  ellipsoidal  coordinates  and  wind-hub  coordinates  are  derived. 

If  we  denote  the  x*w  as  index  representation  of  wind-hub  coordinates  and  x3u  as  the 
index  notation  of  ellipsoidal  coordinates,  we  have  dxlw/dxl  =  dxw/dv,  dxlw/dx2v  = 
dxw/drj,  •  •  •,  dx\ / dx\  —  dyw/di>,  and  etc..  Thus,  we  have 


(3.3) 


This  can  be  further  written  as 
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where  s,  =  ^(1  +  h2)(l  —  L/2)  and  s2  =  v2  +  rf 
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We  also  have  the  relationship  between  the  hub-wind  and  relative  wind  coordinates, 
Fig.  3.1. 
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where  xi  represent  relative  wind  coordinates. 


3.1.2  Gradient  of  Pressure  Function 


The  pressure  function,  Eq.  (2.11),  is  attained  in  the  ellipsoidal  coordinates.  For  in¬ 
duced  velocity  calculation,  we  need  the  gradient  of  the  pressure  function  with  respect 
to  the  relative  wind  frame  and  wind-hub  frame. 

Applying  the  chain  rule,  we  attain  the  gradient  of  the  pressure  function  with 
respect  to  the  relative  wind  frame. 
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(3.6) 


Similarly,  the  gradient  of  the  pressure  function  with  respect  to  the  wind-hub  frame 

can  be  obtained  as  „ 
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(3.7) 


Integration  of  the  pressure  gradient  along  the  free  stream  line  results  in  all  three 
components  of  the  induced  velocity  associated  with  lifting  rotors.  The  results  are 
given  in  the  following  section. 


3.1.3  Three  Components  of  Induced  Velocity 

To  calculate  all  three  components  of  the  induced  velocity  at  an  arbitrary  flowfield 
point,  the  integrations  of  the  spatial  derivatives  of  the  pressure  function  need  to  be 
performed.  Based  on  the  finite  state  dynamic  wake  theory,  these  integrals  are  the 
functions  of  the  geometric  relationship  between  the  rotor  and  the  flowfield  point  of 
interest  for  interference.  The  induced  velocities  thus  obtained  are  both  nonuniform 
and  unsteady.  Moreover,  the  interference  model  can  be  truncated  at  various  frequency 
ranges  in  the  interest  of  different  applications. 

To  make  use  of  the  computational  advantage  of  using  wind-hub  or  relative  wind 
reference  frames,  we  like  to  express  the  induced  velocity  calculations  in  both  of  those 
frames.  In  rotorcraft  analysis  and  simulation  applications,  these  induced  velocities 
are  then,  in  general,  transformed  to  the  nonrotating-hub  frame,  and  further  expressed 
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in  the  local  frame  attached  at  each  aerodynamic  panel.  The  transformations  to  the 
nonrotating  hub  frame  can  be  accomplished  through  the  frame  relationships. 

Component 

The  x-component  of  the  induced  velocity  component  is  calculated  in  the  relative-wind 
frame  due  to  the  availability  of  the  analytical  closed  form.  From  Eq.  (3.1),  we  have 

(3-8) 

Outside  the  rotor  wake,  is  continuous,  and  $v|°°  tends  to  zero.  Therefore 

g*t  (3-9) 

Inside  the  wake,  is  discontinues  wherever  the  integration  crosses  the  rotor 
plane.  Then 

=  —  [^(z^o- ,  y$,  ^cU=o-) 

-  $V(xzk=o+,yt;,zz\i;=o+) 

-  (3.10) 

In  numerical  computation,  a  finite  value  denoted  as  e  needs  to  be  used  for  the 
0”  or  0+.  This  value  is  designated  as  the  “singularity  clearance  distance”  in  the 
numerical  implementation  of  the  theory.  This  number  may  be  seen  as  an  equivalence 
of  the  conventional  vortex  core  size.  In  fact,  it  physically  refers  to  the  discontinuity 
of  the  pressure  across  the  rotor  plane. 

yw- Component 

To  reduce  one  level  of  transformation  from  the  relative  wind  frame  to  the  wind-hub 
frame,  the  y-component  of  induced  velocity  is  calculated  in  the  wind-hub  frame.  From 
Eqs.  (3.1  and  3.5),  we  have 


(3.11) 


Although  the  induced  velocity  component  is  expressed  in  the  wind-hub  frame,  it 
should  be  kept  in  mind  that  the  integration  is  performed  along  the  free  stream  line. 
Substituting  Eq.  (2.11)  into  Eq.  (3.11),  the  ^-component  of  induced  velocity  can 


be  obtained  as 
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where  the  qmc„  and  q™  „  are  the  induced  flow  influence  coefficients  and  are  related 
to  each  specific  distribution  of  the  pressure  function.  They  are 

/■co  f) 

C.  =  J  ^{P™(y)Q";(iri)cos(rni>)\dZ, 

(3.13) 

C«  =  jf  ^[PT<u)Qn^V)stnim.i:)\d( 

(3.14) 

The  derivatives  of  the  associated  Legendre  functions  with  respect  to  the  wind-hub 
coordinates  can  be  obtained  through  chain  rule. 
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zw-Component 


Taking  derivative  of  the  pressure  function  in  the  zw  direction  and  integrating  along 


f,  we  have 


(3.19) 


Substituting  Eq.  (2.11)  into  Eq.  (3.19),  the  z^-component  of  induced  velocity  is 


attained  as 
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where  the  q™„  and  o™  are  the  induced  flow  influence  coefficients  and  can  be  ealeu- 

J-Znj  71  *zw 

lated  from  each  distribution  of  the  pressure  function. 
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As  formulated  in  the  previous  sections,  the  pressure  distribution  coefficients,  ( 
t™cV  and  t™sV  )  are  needed  in  the  calculation  of  the  rotor  induced  velocity  at  arbi¬ 
trary  flowfield  points.  These  pressure  distribution  coefficients  are  obtained  from  the 
rotor  solutions,  Eqs.  2.29  and  2.30 


3.2  Results  and  Discussion 

The  following  sections  examine  the  induced  velocity  computed  from  the  finite  state 
dynamic  wake  interference  model  and  investigates  its  effects  on  helicopter  trim  and 
control  responses. 


3.2.1  Induced  Flowfield 

The  main  rotor  induced  velocity  at  the  empennage  (horizontal  and  vertical  surfaces) 
are  computed  for  the  UH-60  helicopter  at  86  Knots,  Figs.  3.2  and  3.3.  The  results 
are  also  compared  with  the  predictions  of  a  ‘flat’  prescribed  wake  model,  Ref.  [45], 
and  free  wake  model,  Ref.  [15].  Figure  3.2  is  the  predictions  of  the  induced  velocity 
distribution  at  the  location  of  the  horizontal  stabilizer  of  UH-60  helicopter  at  n  =  0.2 
and  Ct  =  0.007.  The  induced  inplane,  sidewash,  and  normal  components  shown  in 
the  rotor  hub  frame  are  normalized  by  the  main  rotor  uniform  downwash  at  the  rotor 
plane.  The  UH-60  horizontal  surface  has  a  full  span  of  about  0.5 R.  The  finite  state 
wake  model  prediction  has  a  close  agreement  with  the  free  wake  model,  Ref.  [15]. 
The  well  established  fact  that  the  downwash  at  the  advancing  side  is  stronger  than  at 
the  retreating  side  is  well  predicted  by  all  the  three  models.  But,  the  flat  wake  model 
overestimates  about  30%  of  the  peak  in  both  the  sidewash  and  normal  components. 
There  is  no  data  available  from  the  free  wake  model,  Ref.  [15],  for  the  inplane 
component  prediction.  Comparing  the  flat  wake  and  finite  state  wake  models,  the 
former  underestimates  the  inplane  component.  But,  for  the  calculation  of  airloads  on 
the  horizontal  stabilator,  only  the  normal  induced  velocity  is  of  significance.  Although 
there  is  no  direct  experimental  data  for  correlation,  the  induced  flow  distribution 
predicted  here,  in  general,  compares  well  with  the  experimental  pattern  as  shown  in 
Ref.  [45]. 

The  induced  velocity  distribution  at  the  vertical  tail  location  is  presented  in  Fig¬ 
ure  3.3.  Of  the  three  velocity  components,  the  sidewash  component  directly  changes 
the  angle  of  attack  of  the  tailplane,  and  thus  has  a  much  larger  effect  on  the  tailplane 
aerodynamics  than  the  other  two  velocity  components.  Only  the  sidewash  data  is 
available  from  Ref.  [15]  for  comparison.  It  is  seen  that  the  finite  state  wake  results 
closely  agree  with  those  of  vortex  wake  models.  There  is  a  slight  shift  at  the  crossover 
where  the  sidewash  changes  its  direction.  This  is  probably  due  to  either  the  incon¬ 
sistency  in  the  vertical  fin  location  between  the  two  calculations  or  the  trim  setting 
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difference. 

One  of  the  distinctive  features  of  the  finite  state  rotor  wake  model  is  its  capability 
to  predict  the  unsteady  and  nonuniform  interference  induced  velocity.  The  interfer¬ 
ence  solution  from  the  finite  state  dynamic  wake  provides  both  the  time  averaged 
induced  flow  effect  and  the  unsteady  time  variation  of  the  induced  velocity  on  the 
impinged  aerodynamic  surfaces.  A  good  example  for  understanding  this  effect  is  to 
consider  a  trailing  vortex  with  time  varying  strength  from  the  main  rotor  that  takes 
time  to  reach  the  aerodynamic  surface  and  then  passes  it.  This  would  create  a  time 
varying  history  of  the  wake  induced  velocity  on  the  aerodynamic  surface.  Thus,  its 
impact  will  affect  both  the  magnitude  and  phase  of  the  helicopter  responses.  The 
modeling  of  this  unsteady  wake  effect  is  a  significant  step  toward  the  improvement  of 
the  fidelity  of  rotorcraft  flight  simulations. 

Figure  3.4  is  the  time-varying  induced  velocity  at  the  fuselage  aerodynamic  center 
and  the  horizontal  tail  location  at  40  Knots  and  80  Knots  forward  flight  speeds.  The 
induced  velocity  varies  at  a  frequency  of  4/rev  corresponding  to  the  number  of  rotor 
blades. 

3.2.2  Effects  on  Helicopter  Trim 

The  effects  of  rotor  interference  on  the  helicopter  trim  are  examined  in  Fig.  3.5 
where  the  trim  results  for  the  UH-60  helicopter  in  level  forward  flight  are  shown.  The 
comparison  of  the  helicopter  trim  is  made  with  and  without  the  main  rotor  wake 
interference,  As  seen,  the  main  rotor  interference  is  most  significant  at  forward  flight 
speeds  around  40  Knots.  It  is  found  that  the  large  change  in  the  longitudinal  stick 
position  and  helicopter  pitch  attitude  is  due  to  the  impingement  of  main  rotor  wake 
on  the  horizontal  tail.  Figure  3.5  also  compares  the  trim  results  of  the  three  state 
inflow  interference  model  with  the  empirical  interference,  Ref.  [46].  A  close  agree¬ 
ment  for  most  control  channels  and  fuselage  attitudes  between  the  two  interference 
models  is  obtained.  A  discrepancy  is  noticed  for  the  longitudinal  control.  More  for¬ 
ward  longitudinal  stick  is  required  for  the  finite  state  wake  interference  than  for  the 
empirical  interference.  This  discrepancy  is  improved  when  a  six  state  dynamic  wake 
interference  is  used. 

3.2.3  Effects  on  Helicopter  Response 

The  UH-60  response  at  60  Knots  to  a  one  inch  aft  longitudinal  stick  is  examined  in 
Fig.  3.6.  The  prediction  of  the  vehicle  response  from  three  models  are  presented  and 
compared  with  flight  test  data.  The  first  model  uses  3-state  induced  inflow  with  its 
interference  on  the  fuselage,  horizontal,  vertical  tails  and  the  tail  rotor.  The  second 
model  uses  a  uniform  inflow  state  and  its  interference.  Between  the  uniform  inflow 
model  and  the  3-state  inflow  model,  there  are  differences  in  both  the  interference 
modeling  and  main  rotor  aerodynamics.  To  examine  the  interference  effect,  the  results 
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from  another  model  that  uses  the  3-state  inflow  model  for  the  main  rotor,  but  allows 
for  interference  from  the  uniform  inflow  state  is  also  presented  for  comparison.  In 
the  last  model,  the  interference  effect  is  related  to  the  rotor  thrust  only.  No  effect  of 
rotor  hub  moments  on  the  interference  is  considered. 

Compared  to  the  uniform  main  rotor  inflow  model,  the  3-state  inflow  model  im¬ 
proves  correlation  in  both  the  on-axis  (pitch  rate)  and  off-axis  responses.  The  inter¬ 
ference  effect  on  the  yaw  response  is  remarkable.  The  3-state  inflow  with  interference 
model  has  the  best  agreement  with  the  flight  test  data,  Ref.  [47]. 

Examining  the  interference  velocity  at  the  tail  rotor,  Fig.  3.7,  reveals  a  large 
variation  in  the  sidewash  induced  flow  component  predicted  by  the  3-state  inflow 
interference  model.  This  variation  can  not  be  present  in  a  table  look-up  method 
due  to  the  dynamic  change  of  relative  position  between  the  main  rotor  wake  and 
interference  location  during  the  maneuvering  flight. 

The  variation  of  this  sidewash  induced  velocity  has  a  significant  impact  on  the 
changes  in  the  angle  of  attack  of  tail  rotor  blade  sections.  The  angle  of  attack  changes 
further  result  in  a  corresponding  variation  in  the  yaw  torque  that  the  tail  rotor  applies 
at  the  helicopter  center  of  gravity,  Fig.  3.7. 


-  40  Knots 

- 80  Knots 


Oscillating  Induced  Velocity  at  Fuselage 


Azimuth  (deg) 


Oscillating  Induced  Velocity  at  Horizontal  Stabilator 


Figure  3.4:  Time  Varying  Induced  Velocity  at  Fuselage  and  Horizontal  Tail 


Chapter  4 

Finite  State  Inflow  in  Vortex  Ring 
State 


As  known,  all  analytical  models  break  down  at  the  vortex  ring  state.  One  of  the  tasks 
under  this  contract  is  to  extend  the  finite  state  wake  model  to  the  vortex  ring  state. 
To  accomplish  this,  two  important  issues  need  to  be  addressed. 


4.1  Vortex  Ring  State  Boundary 


The  first  is  to  define  the  boundary  between  the  vortex  ring  and  the  normal  flow 
states.  A  vortex  ring  state  boundary  criteria  is  required  for  a  smooth  transition  of 
the  finite  state  dynamic  wake  model  from  the  normal  flow  state  to  the  vortex  ring 
condition.  In  practice,  a  validated  analytical  vortex  ring  state  boundary  definition  is 
invaluable  as  a  guidance  for  a  safe  helicopter  flight.  When  a  helicopter  is  descending 
vertically  or  in  a  steep  descent  angle,  a  definite  slipstream  ceases  to  exist.  The  rotor 
downwash  and  normal  free  stream  are  in  the  opposite  flow  directions.  This  causes  the 
rotor  wake  be  pushed  close  to  the  rotor  plane,  and  the  flow  around  the  rotor  becomes 
highly  turbulent.  This  flow  condition  is  usually  called  rotor  vortex  ring  state.  Based 
on  physical  observations,  the  rotor  vortex  ring  state  entrance  condition  occurs  when 
the  value  of  the  flow  component  of  total  relative  air  in  the  negative  rotor  downwash 
direction  is  less  than  a  certain  value,  i.e., 
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where  Vt  is  the  total  relative  air  velocity  at  the  rotor  plane  and  Wio  is  the  induced 
velocity.  The  criteria  value  where  the  rotor  enters  the  vortex  ring  state,  •uCr,  can  be 
further  expressed  as 

Vcr  =  9hWih  (4-2) 
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where  Wih  is  the  rotor  downwash  in  hover  and  g h  can  be  attained  from  the  measured 
data,  Refs.  [23]  to  [26].  In  implementation,  this  criteria  is  built  into  the  coefficients 
of  the  finite  state  dynamic  wake  equations  so  that  a  smooth  transition  to  the  vortex 
ring  state  is  obtained. 


4.2  Vortex  Ring  Model 

Rotor  total  inflow  is  defined  as  the  sum  of  the  rotor  downwash  and  normal  free  stream, 
A  zero  rotor  total  inflow  can  result  from  descent  flight  for  a  main  rotor  or  sideward 
flight  for  a  tail  rotor.  The  zero  total  inflow  implies  a  zero  air  mass  flow  through  the 
rotor  in  descent  flight.  The  existing  finite  state  dynamic  wake  model  has  a  singularity 
at  the  zero  air  mass  flow  condition  and  is  ill-conditioned  near  the  zero  air  mass  flow 
state.  To  overcome  the  singularity  and  improve  the  ill-condition  in  the  vortex  ring 
state,  a  modification  of  the  mass  flow  parameter  in  the  finite  state  wake  equations  is 
required. 

The  vortex  ring  state  is  highly  turbulent.  The  air  viscosity  effect  plays  a  signif¬ 
icant  role  on  the  air  flow  behavior.  No  reliable  analytical  fluid  dynamics  modeling 
is  available  at  this  time.  The  approach  here  is  to  extend  the  finite  state  wake  flow 
parameter  into  the  vortex  ring  state  based  on  a  semi-empirical  model  derived  in  the 
following. 

We  have  experimented  with  a  decay  vortex  length  model,  Refs.  [28]  and  [16],  in 
the  generalization  of  the  uniform  inflow  model.  Since  the  model  is  first  derived  for 
vertical  descent  only,  the  transition  between  normal  flight  and  the  vortex  ring  state 
in  forward  flight  is  not  smooth.  This  model  also  introduces  an  additional  nonlinear 
equation  associated  with  the  decay  vortex  length.  Dependent  upon  the  vortex  ring 
state  entrance  point  there  could  be  a  convergence  problem  in  attaining  the  solution. 
For  this  model  to  work  properly  in  both  axial  and  forward  flight,  a  further  modification 
is  needed. 

There  is  an  alternative  to  generalize  the  finite  state  wake  equations  in  the  vortex 
ring  state.  The  singularity  in  the  flow  parameter  of  the  finite  state  wake  equation 
can  be  removed  by  fitting  the  data  from  the  universal  induced  flow  curve,  Ref.  [27], 
through  the  vortex  ring  region.  To  assure  a  smooth  transition  between  the  universal 
induced  flow  data  and  theoretical  prediction  outside  the  vortex  ring  state,  we  match 
them  along  the  vortex  ring  state  entrance  boundary,  Eq.  (4.1). 

Under  normal  flight  conditions,  the  flow  parameter  is  calculated  from  the  total 
resultant  air  flow  at  the  rotor  plane  as  follows, 

VT  “  '/*T^  (4'3) 

where  /z  is  advance  ratio,  and  A  is  the  total  inflow.  To  derive  the  parameter  applicable 
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in  the  vortex  ring  region,  we  assume  that 

77"  =  hh  (*■*) 

*ring 

The  parameter  J\  is  formulated  for  vertical  descent  and  can  be  written  as 

j\  =  fring[c\{^f  +  Wih)2  +  d\{Xf  +  Wih)  +  e\]  (4.5) 

where  w^  is  the  nondimensional  momentum  theory  value  of  induced  flow  in  hover 
and  A/  is  the  free  stream  air  flow  nondimensionalized  by  rotor  tip  speed.  The  /rins  is 
derived  from  the  universal  induced  flow  data  and  is  approximated  by  a  cubic  function 
as 

fnn,  =  +  Mf-)2  +  Cr(^-)+dr\  (4.6) 

Ct  w^  w^ 

A  curve  fit  of  the  universal  induced  flow  data  in  the  range  of  —2  <  Xf/wih  <  —0.5 
gives  ar  =  0.9606,  bT  =  2.4168,  c?  =  1.2666,  and  dr  =  1.7099. 

For  the  generalized  induced  flow  formulation  in  axial  flight,  Eq  (4.5),  there  are 
three  coefficients  c\,  d\ ,  and  e\  to  be  defined.  They  can  be  determined  by  matching 
the  theoretical  induced  flow  prediction  with  the  universal  data  on  the  the  vortex  ring 
state  boundary.  Three  matching  points  are  selected.  These  include  vortex  ring  state 
entrance  and  exit  points  in  vertical  descent  and  a  central  point  inside  the  vortex  ring 
region.  At  vortex  ring  state  entrance,  we  have 


Xf  =  -Vcr  =  -ghWih 

(4.7) 

Vh»„  =  Vt  =  ^(-gh  +  V^+sib 

(4.8) 

At  vortex  ring  state  exit  to  the  windmill  state,  we  have 

Xf  =  -2 

(4.9) 

Vring  ~  Wih 

(4.10) 

At  vortex  ring  state  central  point, 

Xf  =  -wih 

(4.11) 

V  —  * 

v  Ting  —  r 

Jring 

(4.12) 

Combining  the  above  three  conditions,  the  coefficients  can  be  determined 


—  1  +  (1  —  gft)(A  —  1) 
(1  -  gh)(2  -  gh)w;k 
B  -  1  -  (1  -  gh)2(A  -  1) 

(1  -  s»)( 2  -  g^Wik 


(4.13) 

(4.14) 

(4.15) 
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where 


A 

B 


1 

fring^ih 

2A 


~9h  + 


(4.16) 

(4.17) 


The  parameter  J\  approximates  the  vertical  descent.  The  parameter  is  derived 
for  forward  flight  descent. 


Jfj.  =  —  v)2  +  —  /X)  +  kp  (4-18) 


First,  to  recover  the  axial  flight  modeling,  we  set  to  unity,  resulting  in 

9n  +  hp  +  kft  —  1 


(4.19) 


The  second  condition  is  to  have  VTing  match  VT  for  normal  flight,  Eq.  (4.3),  and  on 
the  transition  boundary,  Eq.  (4.1).  Let  ^  denote  the  criteria  value  of  advance  ratio 
beyond  which  no  vortex  ring  state  can  occur  regardless  of  helicopter  descent  rate. 
The  parameter  ^  is  determined  from  the  vortex  ring  state  criteria  Eq.  (4.1).  This 
condition  yields 

tV  =  \fH  +  (V  +  >"io)2  (4.20) 

where  Wio  is  the  momentum  theory  value  of  induced  flow  at  the  boundary 


C?  1 

2  \jy-2  +  (A/  +  Wut)2 


(4.21) 


The  last  condition  is  to  match  the  universal  empirical  data  at  one  selected  forward 
descent  condition  inside  the  vortex  ring  region.  The  point,  /u  =  has  been 

selected  for  the  fitting  since  it  is  located  midway  between  the  vertical  and  forward 
descent  boundary.  Thus,  we  have 


fring[9n0-  nWih)  +  h^(l  nWih)  +  fJJ  \2  >  ( \  ~7~7~t  Tv; 

2  2  V(5^/,)2  +  (A/  +  fJringWihr 


(4.22) 


where  the  parameter  /M  is  obtained  from  the  universal  induced  flow  data.  Solving  the 
combined  equations  (4.19,  4.20,  and  4.23)  gives 


kp  —  1  (Jn  hp 

_  ~  jO  ~  Mcr(l  —  BjX 

2A^cr(/^cr  2)  /icrt/)j/i(l  ^Wih) 

hp  =  [1  Afj,  +  9^l^cr{l^cr  2)] 


(4.23) 

(4.24) 


(4.25) 
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where 


_  1JJio\n=lieT 
'  ■ 

B'  =  WtZ, 

VTficr  =  \/  ^cr  +  (^/  +  wio)2 


(4.26) 

(4.27) 

(4.28) 


4.3  Results  and  Discussion 

Figure  4.1  presents  the  results  from  the  above  formulation.  The  induced  velocity  is 
normalized  by  the  momentum  theory  value  in  hover,  Wih,  and  plotted  as  a  function 
of  vertical  climb/descent  rate.  The  climb/descent  rate  is  also  normalized  by  wih. 
The  upper  figure  is  for  axial  flight  and  the  lower  is  the  result  at  /j.  =  0.5^.  The 
results  are  also  compared  to  the  data  obtained  from  the  universal  inflow  curve,  Ref. 
[27].  As  shown,  the  singularity  in  the  theoretical  formulation  has  been  removed  and 
the  prediction  reasonably  matches  the  measured  data.  The  remaining  discrepancy 
between  the  analytical  formulation  and  data  is  not  a  concern  since  the  agreement  is 
well  within  the  scatter  of  the  measurement. 

Figure  4.2  shows  the  induced  flow  variations  with  both  vertical  and  forward  de¬ 
scent.  A  smooth  transition  between  the  theoretical  and  semi-empirical  formulations 
as  derived  above  has  been  obtained.  With  this  new  formulation,  the  enhanced  fi¬ 
nite  state  induced  flow  model  can  now  be  applied  in  all  flight  ranges  for  rotorcraft 
simulation. 


Normalized  climb  rate 


Figure  4.1:  Comparison  of  induced  flow  between  prediction  and  measured  data 


Chapter  5 
Conclusions 


For  rotorcraft  simulation  model  enhancements,  we  have  developed  a  finite  state 
dynamic  wake  formulation  that  includes  the  wake  distortion  effect  due  to  the  hub 
rotation  to  overcome  the  poor  prediction  of  rotorcraft  off-axis  response.  The  formu¬ 
lation  also  includes  the  rotor  aerodynamic  interference  on  other  lifting  surfaces,  such 
as  the  horizontal  tail  and  the  vertical  fin,  and  the  tail  rotor.  An  enhancement  of 
the  finite  state  dynamic  wake  model  in  the  vortex  ring  state  has  also  been  made  to 
generate  a  simulation  model  applicable  across  the  full  range  of  flight  conditions. 

The  following  conclusions  can  be  made  from  the  Phase  One  effort: 

•  A  significant  breakthrough  has  been  made  in  the  development  of  a  parametric 
formulation  using  wake  distortion  influence  coefficients  due  to  hub  rotation  in 
both  hover  and  forward  flight.  It  is  found  that  the  perturbation  inflow  states 
due  to  the  wake  distortion  vary  linearly  with  the  tip-path-plane  rotation  rate 
during  maneuvering  flight.  They  are  also  only  functions  of  flight  speed  and  tip- 
path-plane  angle  of  attack  for  pitch,  and  flight  speed  and  tip-path-plane  lateral 
tilt  angle  for  roll. 

•  Inclusion  of  a  wake  distortion  effect  due  to  the  hub  rotation  in  the  finite  state  in¬ 
duced  flow  model  has  greatly  improved  the  correlation  of  simulation  results  with 
the  flight  test  data  in  pitch/roll  coupling.  The  improvement  of  the  correlation 
has  been  obtained  in  both  hover  and  forward  flight. 

•  The  wake  distortion  effect  due  to  the  hub  rotation  is  more  pronounced  in  hover 
and  low  speed.  The  hub  rotation  effect  is  mainly  on  the  on-axis  inflow  distribu¬ 
tion  (i.e.,  the  fore-to-aft  inflow  variation  due  to  the  pitch  and  side-to-side  inflow 
change  due  the  roll).  The  cross-axis  inflow  variation  due  to  the  hub  rotation 
(i.e.,  the  side-to-side  inflow  distribution  due  the  pitch,  or  the  fore-to-aft  inflow 
distribution  due  to  the  roll)  is  much  less  significant. 

•  The  capability  of  the  finite  state  dynamic  wake  model  has  been  extended  in 
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modeling  the  aerodynamic  interference  between  the  main  rotor  and  the  empen¬ 
nage.  The  results  have  shown  that  the  aerodynamic  interference  is  an  another 
factor  that  has  an  impact  on  the  cross-axis  coupling  response,  such  as  the  yaw 
due  to  collective  or  longitudinal  inputs,  or  the  pitch  due  to  sideslip,  etc.. 

•  The  enhanced  model  formulation  has  removed  the  singularity  of  inflow  variation 
in  the  vortex  ring  state.  A  smooth  induced  flow  variation  has  been  obtained 
in  both  vertical  and  forward  flight  descent.  The  analysis  results  well  match 
measured  universal  induced  flow  curves. 

•  The  formulation  of  the  vortex  ring  state  model  has  provided  a  criteria  that 
describes  the  vortex  ring  state  boundary.  This  can  be  used  as  a  flight  safety 
guide. 


Chapter  6 
Future  Plan 


We  have  successfully  achieved  the  Phase  One  objectives  for  rotorcraft  simulation 
model  enhancement.  The  Phase  One  results  have  shown  that  the  approach  is  effective 
and  feasible.  The  formulation  developed  can  be  fully  implemented  in  a  comprehensive 
rotorcraft  simulation  program  for  both  real-time  flight  simulation  or  non-real  time 
analysis. 

Further  Phase  Two  model  enhancement  efforts  will  be  made  in  the  following  ma¬ 
jor  areas. 

•  Wake  distortion  modeling:  We  will  fully  implement  the  wake  distortion 
formulation  as  developed  in  Phase  One  in  a  comprehensive  rotorcraft  simulation 
program  to  address  rotorcraft  response  in  both  trim  and  maneuvering  flight. 
Phase  One  has  only  performed  the  feasibility  study.  For  the  formulation  to  be 
fully  validated,  more  cases  will  be  tested  and  correlated  with  measured  data. 

•  Aerodynamic  interaction  modeling:  The  wake  distortion  due  to  the  hub 
rotation  mainly  affects  the  helicopter  pitch/roll  coupling.  In  forward  flight, 
there  are  other  physical  phenomena  that  affect  the  simulation  correlation  with 
the  flight  test  data  in  helicopter  coupling  response.  One  of  them  is  the  fuselage 
induced  fore-to-aft  inflow  variation  at  the  rotor  plane  that  can  influence  the 
helicopter  roll  response  in  forward  flight.  In  Phase  Two,  a  panel  module  for 
interference  velocity  calculation  will  be  developed  and  implemented  to  model 
the  fuselage  and  other  lifting  surfaces  with  low  aspect  ratio. 

•  Special  operation  modeling:  We  will  address  the  simulation  of  the  vortex 
ring  state  in  descent  and  ground  effect  modeling  during  ship  deck  landing.  For 
the  vortex  ring  state  modeling,  the  formulation  developed  in  Phase  One  will  be 
fully  implemented  in  a  comprehensive  simulation  program,  such  as  FLIGHT- 
LAB.  The  vortex  ring  state  boundary  will  also  be  computed  in  the  simulation 
as  a  flight  safety  guide  and  a  set  of  critical  azimuth  control  margin  plots  will 
be  generated  for  a  selected  vehicle.  For  the  ship  deck  landing,  both  dynamic 
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ground  effect  when  the  ship  deck  is  moving  and  partial  ground  effect  will  be 
modeled. 

•  Swirl  velocity:  We  will  enhance  the  finite  state  dynamic  wake  model  to  cal¬ 
culate  the  induced  flow  in  the  azimuthal  direction  due  to  the  in-plane  loading 
known  as  swirl  effect.  The  swirl  is  a  unique  feature  that  is  especially  important 
for  prop-rotors  as  used  in  tilt  rotor  configurations.  Due  to  the  high  blade  pitch 
setting,  the  swirl  velocity  becomes  significant  in  affecting  both  performance  and 
response  analysis  of  the  tilt  rotor. 

•  Fuselage  airloads  calculation:  Most  rotorcraft  simulations  still  rely  on  em¬ 
pirical  data  for  fuselage  airload  calculation.  As  a  first  step,  a  panel  model  will 
be  developed  for  the  computation  of  the  steady  fuselage  airloads.  The  flow 
separation  will  be  considered  by  a  boundary  layer  correction  with  the  potential 
flow  solution.  The  panel  airload  model  will  have  applications  beyond  rotorcraft, 
such  as  automobile  simulation. 
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